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Abstract 

We review several aspects of flavour-diagonal CP violation, focussing on the role 
played by the electric dipole moments (EDMs) of leptons, nucleons, atoms and 
molecules, which consitute the source of several stringent constraints on new CP- 
violating physics. We dwell specifically on the calculational aspects of applying the 
hadronic EDM constraints, reviewing in detail the application of QCD sum-rules 
to the calculation of nucleon EDMs and CP-odd pion-nucleon couplings. We also 
consider the current status of EDMs in the Standard Model, and on the ensuing con- 
straints on the underlying sources of CP-violation in physics beyond the Standard 
Model, focussing on weak-scale supersymmetry. 



1 Introduction 



The search for violations of fundamental symmetries has played a central role 
in the development of particle physics in the 20 th century. In particular, tests 
of the discrete symmetries, charge conjugation C, parity P, and time-reversal 
T, have been of paramount importance in establishing the structure of the 
Standard Model (SM). Perhaps the most famous example was the discovery 
of parity violation in the weak interactions [1], which led to the realization 
that matter fields should be combined into asymmetric left- and right-handed 
chiral multiplets, one of the cornerstones of the Standard Model. The observa- 
tion of CP violation, via the mixing of Kaons [2], also subsequently provided 
strong evidence for the presence of three quark and lepton generations, via 
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the Kobayashi-Maskawa mechanism [3] , prior to direct experimental evidence 
for the third family. 

It is interesting to recall that one of the first tests of this kind, actually pre- 
dating the discovery of parity violation in the weak interactions, was a probe 
of parity invariance within the - at the time unknown - theory of the strong 
interactions. In 1949, Purcell and Ramsey argued, in a way that at the time 
was not fully appreciated that, lacking a theory of the strong interactions, 
there was no way to "derive" parity invariance and thus one must confirm 
parity conservation by experimental tests, or discover the lack thereof, rather 
than rely on a belief that nature "must be symmetric" . As a probe of parity 
violation, Purcell and Ramsey proposed that one consider an intrinsic elec- 
tric dipole moment of the neutron. Placed in a magnetic and electric field, 
a neutral nonrelativistic particle of spin S can be described by the following 
Hamiltonian, 

= -^B.|-dE.| . (1.1) 

Under the reflection of spatial coordinates, P(B-S) = BS, whereas P(E-S) = 
— E • S. The presence of a non-zero d would therefore signify the existence of 
parity violation. It was soon realized that d also breaks time-reversal invari- 
ance. Indeed, under time reflection, T(B • S) = B • S and T(E • S) = — E • S. 
Therefore a non-zero d may exist if and only if both parity and time reveral 
invariance are broken. Analysis of the existing experimental data on neutron 
scattering from spin zero nuclei led to the conclusion, \d n \ < 3 x lCT 18 ecm 
[4]. Such a result probes physics at distances much shorter than the typical 
scale of nuclear froces ~ lfm, or the Compton wavelength of the neutron. This 
initial limit on the neutron EDM implied that P and T were good symmetries 
of the strong interactions at percent-level precision. 

It was only some 25 years later with the emergence of QCD that the possibility 
of T- violation (or CP- violation, on assuming the CPT theorem) in the strong 
interactions had some theoretical underpinning. Indeed, QCD allows for the 
addition of a dimension-four term, known as the #-term, with a dimensionless 
coefficient 9 which, if nonzero, would signify the violation of both P and T. 
This term is somewhat unusual, being a purely topological boundary term, 
but its value determines a superselection sector in QCD [6] and its presence 
is intrinsically tied to an elegant feature of the theory, namely the mechanism 
via which the mass of the rf meson is lifted well-above the scale one might 
naturally expect given its apparent status as a Goldstone boson [5]. However, 
were 9 ~ 0(1), one would predict a neutron EDM of sufficient size to ensure 
that the original analysis of Purcell and Ramsey would have detected it. In fact 
9 is now known to be tuned to better than one part in 10 9 ! This tuning is the 
well-known strong CP problem of the Standard Model, which has been with 
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us for more than 25 years, and has led to interesting dynamical mechanisms 
for its resolution; some of these have important consequences and predictions 
for other aspects of particle physics and cosmology. 

The required tuning of this CP-odd parameter in QCD comes into sharp focus 
when we put QCD into its rightful place within the Standard Model, which 
necessarily means coupling it to the electroweak sector and massive quarks in 
particular. In this case the physical value of 9 acquires a contribution from the 
overall phase of the quark mass matrix. In this sense the strong CP problem 
can be phrased as the absence, to high precision, of flavour- diagonal CP- 
violation within the Standard Model. This situation could not contrast more 
strongly with the situation in the flavour- changing sector, which is where all 
currently observed CP- violating effects reside. Indeed, the original discovery of 
CP- violation in the system of neutral Kaons [2], can be explained within this 
sector through the elegant and indeed rather minimal model of Kobayashi and 
Maskawa, which links CP- violation to the single physical phase in the unitary 
CKM mixing matrix describing transitions between the three generations of 
quarks [3]. This picture has recently received significant support - indeed 
essential confirmation - through experiments using neutral B mesons [7]. In 
contrast to 9, the phase in the CKM mixing matrix requires no tuning at 
all - its effects are nicely masked in the appropriate channels by the flavour 
structure of the Standard Model. Indeed, it turns out that the predictions for 
any CP- violating effect in the flavour-conserving channel induced by CKM 
mixing are minuscule, thus denying any hopes of detecting the experimental 
manifestation of CKM physics in these channels in the foreseeable future. 

Searches for flavour-diagonal CP- violation, while insensitive to the CKM phase, 
thus inherit on the flip-side the status as one of the unique, essentially "back- 
ground" free, probes of new physics. Electric dipole moments, through contin- 
uous experimental development since the work of Purcell and Ramsey, remain 
our most sensitive probes of this sector. All existing searches have failed to 
detect any intrinsic EDM, and indeed the precision to which EDMs are now 
known to vanish is remarkable, and sufficient to render them some of the 
most important precision tests of the Standard Model. In this more general 
context, the strong CP problem, associated with the tuning of 9, becomes 
just the most highly tuned example among many possible CP-odd operators 
which could contribute to the observable EDMs of nucleons, leptons, atoms 
and molecules. Anticipating the presence of such CP-odd sources is not with- 
out motivation. Indeed, one of the strongest motivations comes from cosmol- 
ogy, where the success of the inflationary scenaro, together with the observed 
cosmological dominance of baryons over antibaryons, suggests that a non-zero 
baryon number was generated dynamically in the early Universe. According 
to the Sakharov criteria [8], this requires a source of CP- violation, and within 
the Standard Model, neither the Kobayashi-Maskawa phase nor the #-term 
can create conditions that would lead to the generation of an appreciable net 
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Table 1 



Current constraints within three representatve classes of EDMs 



Class 


EDM 


Current Bound 


Paramagnetic 
Diamagnetic 
Nucleon 


205 T/ 
199 Hg 

n 


\d Tl \ < 9 x 10- 25 ecm (90% C.L.) [10] 
\d Ug \ < 2 x 10- 28 ecm (95% C.L.) [11] 
\d n \ < 6 x 10- 26 ecm (90% C.L.) [9] 



baryon number. This strongly suggests the presence of another, yet to be dis- 
covered, source of CP violation in nature. Although exceptions exist, e.g. the 
leptogenesis scenario, EDMs generally provide a highly sensitive diagnostic for 
these new CP-odd sources. 

The second prominent motivation arises from theoretical prejudices about 
the physics of the Fermi scale, i.e. the mechanism for electroweak symmetry 
breaking, currently the focus of intense theoretical and experimental work. 
There are several theoretical motivations to believe that new physics, beyond 
the SM Higgs boson, should become apparent at, or just above, this scale, 
with weak-scale supersymmetry (SUSY) being a prominent example. Flavour- 
diagonal CP- violation constitutes a powerful probe of these scales, since any 
new physics need not provide the same flavour-dependent suppression factors 
as does the SM, while the SM itself constitutes a negligible background. These 
precision tests are thus highly complementary to direct searches at colliders. A 
rough estimate, based on the decoupling of new physics as the inverse square of 
its characteristic energy scale A, currently gives us the possibility to probe an 
order one CP-violating source at up to A ~ 10 6 GeV. In many weakly coupled 
theories, such as SUSY, this scale is somewhat lower, but often is still beyond 
the reach of existing and/or projected colliders. As with the link between the 
Kobayashi-Maskawa mechanism and the three-generation structure, one might 
hope that flavour- diagonal CP- violation, or perhaps the lack thereof, will tell 
us something profound about the matter sector. 

The level of experimental precision achieved in EDM searches has improved 
dramatically since the early work of Purcell and Ramsey, and has been broad- 
ened to many atomic and nuclear quantities. Indeed, following significant 
progress throughout the past decade, the EDMs of the neutron [9], and of 
several heavy atoms and molecules [10-15] have been measured to vanish to 
remarkably high precision. From the present standpoint, it is convenient to 
classify the EDM searches into three main categories, distinguished by the 
dominant physics which would induce the EDM, at least within a generic 
class of models. These categories are: the EDMs of paramagnetic atoms and 
molecules; the EDMs of diamagnetic atoms; and the EDMs of hadrons, and 
nucleons in particular. For these three categories, the experiments that cur- 
rently champion the best bounds on CP- violating parameters are the atomic 
EDMs of thallium and mercury and that of the neutron, as listed in Table 1. 
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The upper limits on EDMs obtained in these experiments can be translated 
into tight constraints on the CP- violating physics at and above the electroweak 
scale, with each category of EDM primarily sensitive to different CP-odd 
sources. For example, the neutron EDM can be induced by CP violation in 
the quark sector, while paramagnetic EDMs generally result from CP violat- 
ing sources that induce the electron EDM. Despite the apparent difference in 
the actual numbers in Table 1, all three limits on d n , d^i, and du g actually 
have comparable sensitivity to fundamental CP violation, e.g. superpartner 
masses and CP- violating phases, and thus play complementary roles in con- 
straining fundamental CP-odd soources. This fact can be explained by the 
way the so-called Schiff screening theorem [16] is violated in paramagnetic 
and diamagnetic atoms. The Schiff theorem essentially amounts to the state- 
ment that, in the nonrelativistic limit and treating the nucleus as pointlike, the 
atomic EDMs will vanish due to screening of the applied electric field within a 
neutral atom. The paramagnetic and diamagnetic EDMs result from violations 
of this theorem due respectively to relativistic and finite-size effects, and in 
heavy atoms such violation is maximized. For heavy paramagnetic atoms, i.e. 
atoms with non-zero electron angular momentum, relativistic effects actually 
result in a net enhancement of the atomic EDM over the electron EDM. For 
diamagnetic species, the Schiff screening is violated due to the finite size of 
the nucleus, but this is a weaker effect and the induced EDM of the atom is 
suppressed relative to the EDM of the nucleus itself. These factors equilibrate 
the sensitivities of the various experimental constraints in Table 1 to more 
fundamental sources of CP violation. 

In this paper we will review in detail the calculational aspects of applying the 
current bounds on EDMs to constrain new CP-violating sources. In order to 
make the discussion as systematic as possible, we will proceed by working our 
way upwards in energy scale, using several effective CP-odd Lagrangians at 
the relevant thresholds. In the next section, we begin by discussing the current 
status of the experimental constraints within three generic classes, namely the 
neutron EDM, and the EDMs of paramagnetic and diamagnetic atoms, and 
describing the contributions to these EDMs at the nuclear scale. We then move 
to the QCD scale and introduce an effective CP-odd effective quark-gluon La- 
grangian which plays an important role in the subsequent analysis. The leading 
term in this effective theory, of dimension four, is the #-term and we breifly 
review the strong CP-problem and some of its proposed resolutions. We then 
turn in Section 3 to QCD computations of the EDMs, and dwell on some of 
the calculational aspects which are currently some of the major sources of 
uncertainty in the application of EDM constraints (for a detailed discussion 
of many other aspects we refer the reader to [17]). In Section 4, we turn to the 
generation of these observables within specific models of CP-violation, review- 
ing first the significant sources of suppression within the Standard Model, and 
then focussing on weak scale supersymmetry, and the MSSM in particular, 
as the source of new physics at the electroweak scale. We discuss the generic 
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constraints that EDMs impose on combinations of CP-violating parameters 
in the SUSY-breaking sector, and also explore some additional effects which 
may arise in special parameter regimes. We also emphasize the stringent EDM 
constraints on combined sources of CP- and flavour- violation in more general 
models. Finally, we conclude in Section 5 with an outlook on future experi- 
mental and theoretical developments. 



2 EDMs as probes of CP violation 

The majority of EDM experiments are performed with matter as opposed to 
anti- matter. Therefore, the conclusion about the relation between d and CP 
violation relies on the validity of the CPT theorem. The interaction c?E • S for 
a spin 1/2 particle then has the following relativistic generalization 

#T,p_odd = ~dE • | — C = -d-^^F^. (2.2) 

Parenthetically, it is worth remarking that the precision of EDM experiments 
has now reached a level sufficient to provide competetive tests of CPT invari- 
ance, since one can also consider a CP-even, but CPT-odd, relativisitic form 
of dE • S, namely C = dt/j^^ipF^n" , with a preferred frame n u = (1, 0, 0, 0), 
which spontaneously breaks Lorentz invariance and CPT. 

The problem of calculating an observable EDM from the underlying CP vi- 
olation in a given particle physics model can be conveniently separated into 
different stages, depending on the characteristic energy/momentum scales. At 
each step the result can be expressed as an effective Lagrangian in terms of 
the light degrees of freedom with Wilson coefficients that encode information 
about CP violation at higher energy scales. As usual in effective field theory, it 
is convenient to classify all possible effective CP violating operators in terms 
of their dimension, with the operators of lowest dimension usually leading to 
the largest contributions. This logic may need to be refined if symmetry re- 
quirements imply that certain operators are effectively of higher dimension 
than naive counting would suggest. This is actually the case for certain EDM 
operators due to gauge invariance, as discussed in more detail below. 

We will present this analysis systematically in order of increasing energy scale, 
working our way upwards in the dependency tree outlined in Fig. 1, which 
allows us to remain entirely model-independent until the final step where 
some high-scale model of CP violation can be imposed and then subjected to 
EDM constraints. 
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Fig. 1. A schematic plot of the hierarchy of scales between the CP-odd sources 
and three generic classes of observable EDMs. The dashed lines indicate generically 
weaker dependencies. 



2.1 Observable EDMs 



Let us begin by reviewing the lowest level in this construction, namely the 
precise relations between observable EDMs and the relevant CP-odd operators 
at the nuclear scale. At leading order, such effects may be quantified in terms 
of EDMs of the constituent nucleons, d n and d p (where the neutron EDM 
is already an observable), the EDM of the electron d e , and CP-odd electron- 
nucleon and nucleon-nucleon interactions. In the relevant channels these latter 
interactions are dominated by pion exchange, and thus we must also consider 
the CP-odd pion-nucleon couplings g W NN which can be induced by CP-odd 
interactions between quarks and gluons. To be more explicit, we write down 
the relevant CP-odd terms at the nuclear scale, 

^nuclear _ £ edm _j_ C^NN + C-eN, (2.3) 

which can be split into terms for the nucleon (and electron) EDMs, 



C edm = - % - E diMFa)^, (2.4) 

i=e,p,n 

the CP-odd pion nucleon intercations, 

C*nn =g ( ^ N Nr a Nn a + £l N NNn° 

+£l N (Nr a N7r a -3Nt 3 Nti°), (2.5) 
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and finally CP-odd electron-nucleon couplings, 



C eN = cfei^eNN + CfeeNi^N + C^'e^ea^eNa^N 

+C { s ) ei-f 5 eNr 3 N + C^eeNi^N + C^e^ea^eNa^^N. (2.6) 

In certain rare cases, CP-odd nucleon-nucleon forces are not mediated by 
pions, in which case the effective Lagrangian must be extended by a variety 
of contact terms e.g. NNNi^y^N, and the like. 

The dependence of the observable EDMs on the corresponding Wilson coef- 
ficients relies on atomic and nuclear many-body calculations which would go 
beyond the scope of this review to cover here (see the reviews [17,18] for fur- 
ther details). However, we will briefly summarize the current status of these 
calculations, before turning to our major focus which is the calculation of these 
coefficients in terms of higher scale CP-odd sources. 

As alluded to earlier on, it is convenient to split the discussion into three 
parts, corresponding roughly to the three classes of observable EDMs which 
currently provide constraints at a similar level of precision; namely: EDMs 
of paramagnetic atoms and molecules, EDMs of diamagnetic atoms, and the 
neutron EDM. 



• EDMs of paramagnetic atoms - thallium EDM 



Paramagnetic systems, namely those with one unpaired electron, are primarily 
sensitive to the EDM of this electron. At the nonrelativistic level, this is far 
from obvious due to the Schiff shielding theorem which implies, since the atom 
is neutral, that any applied electric field will be shielded and so an EDM of the 
unpaired electron will not induce an atomic EDM. Fortunately, this theorem 
is violated by relativistic effects. In fact, it is violated strongly for atoms with 
a large atomic number, and even more strongly in molecules which can be 
polarised by the applied field. For atoms, the parameteric enhancement of the 
electron EDM is given by [19,20,18], 



Z 3 a 2 

rfpara(4) ~ 10 j(j + i/ 2 )(j + r/ e ' (2 - 7) 

up to numerical 0(1) factors, with J the angular momentum and Z the atomic 
number. This enhancement is significant, and for large Z, the applied field 
can be enhanced be a factor of a few hundred within the atom. This feature 
explains why atomic systems provide such a powerful probe of the electron 
EDM, since the "effective" electric field can be much larger than one could 
actually produce in the lab. 
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Although the electron EDM is the predominant contributor to any paramag- 
netic EDM is most models, one should bear in mind that other contributions 
may also be significant in certain regimes. In particular, significant CP-odd 
electron-nucleon couplings may also be generated, due for example to CP vi- 
olation in the Higgs sector. Among these couplings, Cs plays by far the most 
important role for paramagnetic EDMs because it couples to the spin of the 
electron and is enhanced by the large nucleon number in heavy atoms. 

Among various paramagnetic systems, the EDM of the thallium atom cur- 
rently provides the best constraints on fundamental CP violation. A number 
of atomic calculations [20-22] (see also Ref. [17] for a more complete list) have 
established the relation between the EDM of thallium, d e , and the CP-odd 
electron-nucleon interactions Cs'- 

d T1 = -5854 - e 43 GeV x (C^ - 0.2Cj } ). (2.8) 

with Cs expressed in isospin components. The relevant atomic matrix elements 
are known to within 10 — 20% [18]. 

As we discuss later on, current experimental work is focussing on the use of 
paramagnetic molecules, e.g. YbF and PbO [15,23], which can provide an even 
larger enhancement of the applied field due to polarization effects, have better 
systematics, and may bring significant progress in measuring/constraining d e 
and Cs- 



• EDMs of diamagnetic atoms - mercury EDM 



EDMs of diamagnetic atoms, i.e. atoms with total electron angular momentum 
equal to zero, also provide an important test of CP violation [17]. In such 
systems the Schiff shielding argument again holds to leading order. However, 
in this case it is violated not by relativistic effects but by finite size effects, 
namely a net misalignment between the distribution of charge and EDM (i.e. 
first and second moments) in the nucleus of a large atom (see e.g. [18] for a 
review). However, in contrast to the paramagnetic case, this is a rather subtle 
effect and the induced atomic EDM is considerably suppressed relative to the 
underlying EDM of the nucleus. 

To leading order in an expansion around the point-like approximation for 
the nucleus, the contributions arise from an octopole moment (which is only 
relevant for states with large spin, and will not be relevant for the cases con- 
sidered here), and the Schiff moment S, which contributes to the electrostatic 
potential, 
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V E = AtvS- f5(f). 



(2.9) 



CP-odd nuclear moments, such as S, can arise from intrinsic EDMs of the 
constituent nucleons and also CP-odd nucleon interactions. It turns out that 
the latter source tends to dominate in diamagnetic atoms and thus, since 
such interactions are predominantly due to pion exchange, we can ascribe the 
leading contribution to CP-odd pion nucleon couplings g^ NN for i = 0,1,2 
corresponding to the isospin. 

There are of course various additional contributions, which are generically 
subleading, but may become important in certain models. Schematically, we 
can represent the EDM in the form 

ofaia = ddin(S[g n NN, d>N\, C5, Cp, Ct, d e ), (2-10) 

where we note that electro n-nucleon interactions may also be significant, as is 
the electron EDM itself [17] (although in practice the electron EDM tends to 
be more strongly constrained by limits from paramagentic systems and thus 
is often neglected). 

Currently, the strongest constraint in the diamagnetic sector comes from the 
bound on the EDM of mercury - at the atomic level, this is in fact the most 
precise EDM bound in existence. As should be apparent from the above discus- 
sion, computing the dependence of c?Hg o n the underlying CP-odd sources is a 
nontrivial problem requiring input from QCD and nuclear and atomic physics. 
In particular, the computation of S^nn) is a nontrivial nuclear many-body 
problem, and has recently been reanalyzed with the result [24], 

S( 199 Hg) = -0.0004^ (0) - 0.055^ (1) + 0.009^ (2) e fm 3 , (2.11) 

where g = g-^NN is the CP-even pion- nucleon coupling, and = g^ NN denote 
the CP-odd couplings. The isoscalar and isotensor couplings have been sig- 
nificantly reduced relative to previous estimates, while the isovector coupling 
- which generically turns out to be most important - has been less affected 
(within a factor 2). This nonethless provides some indication of the difficul- 
ties inherent in the calculation, and makes precision estimates more difficult. 
Moreover, it is worth noting that the suppression of the overall coefficient in 
front of ggW below 0(0.01) is the result of mutual cancellation between sev- 
eral contributions of comparable size, and therefore is in some sense accidental 
and may not hold in future refinements of these nuclear calculations. 

Putting the pieces together, we can write the mercury EDM in the form, 
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(2.12) 



where we have limited attention to the isovector pion-nucleon coupling and C$ 
which turns out to be the most important for CP violation in supersymmetric 
models. 

• Neutron EDM 

The final class to consider is that of the neutron itself, whose EDM can be 
searched for directly with ultracold neutron technology, and currently provides 
one of the strongest constraints on new CP-violating physics. In this case, 
there is clearly no additonal atomic or nuclear physics to deal with, and we 
must turn directly to the next level in energy scale, namely the use of QCD 
to compute the dependence of d n on CP-odd sources at the quark-gluon level. 
This statement also applies to many of the other quantities we have introduced 
thus far, including in particular the CP-odd pion-nucleon coupling. Indeed, 
it is only paramagnetic systems that are partially immune to QCD effects, 
although even there we have noted the possible relevance of electron-nucleon 
interactions. 

2.2 The structure of the low energy Lagrangian at 1 GeV 

The effective CP-odd flavour-diagonal Lagrangian normalized at 1 GeV, which 
is taken to be the lowest perturbative quark/gluon scale, plays a special role 
in EDM calculations. At this scale, all particles other than the u, d and s 
quark fields, gluons, photons, muons and electrons can be considered heavy, 
and thus integrated out. As a result, one can construct an effective Lagrangian 
by listing all possible CP-odd operators in order of increasing dimension, 

£cff = A}im=4 + Aiim=5 + A}im=6 + ' ' ' • (2-13) 

There is only one operator at dimension 4, the QCD theta term, 



where on account of the axial U(l) anomaly, the physical value of 9 - denoted 
6 - also includes the overall phase of the quark mass matrix, 



dim=4 




(2.14) 



6 = 6 + ArgDetMg. 



(2.15) 
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The anomaly can be used to shuffle contributions between the #-term and 
imaginary quark masses, but only the combination 9 is physical and we choose 
to place it in front of GG taking DetM g to be real. It should be apparent that 
if any of the quarks were massless, we could then rotate 6 away and it would 
have no physical consequences. 

At the dimension five level, there are (naively) several operators: EDMs of 
light quarks and leptons and color electric dipole moments of the light quarks, 



£ d im=5 = Z di ^i( Fa )l^i ~2^ di ^i9s{G^)lb^i, (2.16) 

i=u,d,s,e,fj, i=u,d,s 

where (Fa) and (Go) are a shorthand notation for F^a^ and G a t a a IJ ' v ' . 

In fact, in most models these operators are really dimension-six operators in 
disguise. The reason is that, if we proceed in energy above the electroweak scale 
and assume the system restores SU(2)xC/(l) as in the Standard Model, gauge 
invariance ensures that these operators must include a Higgs field insertion 
[25]. Indeed, were we to write the basis of down quark EDMs and CEDMs 
above the electroweak scale, we should specify the following list of dimension 
six operators [25], 



= ^=Q L [2d™(Ba) + dfV(WV) 

+ df N X a (G a a)] (Q/v)D R + h.c, (2.17) 

which are defined in terms of left-handed doublets Ql = (U, D)l and right- 
handed singlets Dr and the Higgs doublet $, and in terms of the U(l), SU(2), 
and SU(3) field strengths B^, and G^ v . 

The lesson we draw from (2.17) with regard to EDMs is that, if generated, 
these operators must be proportional to the Higgs v.e.v. below the electroweak 
scale, and consequently must scale at least as 1/M 2 for M 3> M w . In practice, 
this feature can also be understood in most models by going to a chiral basis, 
where we see that these operators connect left- and right-handed fermions, 
and thus require a chirality flip. This is usually supplied by an insertion of 
the fermion mass, i.e. df ~ rrif/M 2 , again implying that the operators are 
effectively of dimension six. 

Consequently, for consistency we should also proceed at least to dimension six 
where we encounter the CP-odd three-gluon Weinberg operator and a host of 
possible four-fermion interactions, (ipiTipijiipjiT^^pj), where V denotes several 
possible scalar or tensor Lorentz structures and/or gauge structures, which are 
contracted between the two bilinears. We limit our attention to a small subset 
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of the latter that will be relevant later on, 



£ dim=6 = -w f abc G; u G^ b Gp »> c + £ C y (AAMjilM + ■■■ (2.18) 



In this formula, the operators with CV, are summed over all light fermions. 
Going once again to a chiral basis, we can argue as above that the four-fermion 
operators, which require two chirality flips, are in most models effectively of 
dimension eight. Nonetheless, in certain cases they may be non- negligible. 

2.3 The strong CP problem 

The leading dimension-four term in the CP-odd Lagrangian given in Eq. (2.14) 
has a special status, in that it is a marginal operator, unsuppressed by any 
heavy scale. It is also a total derivative - we can write GG = d^K,^ with /C M 
the Chern-Simons current - and thus plays no role in perturbation theory. 
However, /C M is not invariant under so-called large gauge transformations and 
thus one may expect that the #-term becomes relevant at the nonperturbative 
level. That it does so can be argued at the semi-classical level using instanton 
methods, and more generally can be understood within QCD via this relation 
to the U(l) problem. In particular, we note that the same operator GG arises 
as the #-term in the Lagrangian, and also as an anomaly for the axial U(l) 
current J^, i.e. for massless quarks, 



This leads to an intrinsic link between two physical phenomena: namely the 9- 
dependence of physical quantities, and the absence of a light pseudo-Goldstone 
boson associated with spontaneous breaking of the axial current J% [5] (the 
corresponding state, the rj' is instead rather heavy, m v i ^> m n ). Although it 
would take us too far afield to review the story of this link in detail (see e.g. 
[5,26-30]), let us note that in the large N limit, as discussed by Witten and 
Veneziano [27,28], use of the anomaly equation leads to a simple relation that 
exemplifies this connection, 



where Nf is the number of flavours. This relation expresses the rj' mass in 
terms of the ^-dependence of the vacuum energy in a theory with no light 
quarks. 



W = t^G-G' 



(2.19) 




(2.20) 
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In turn, if we now take for granted that m v > ^> m^, use of the anomaly relation 
allows precise calculations of the ^-dependence of physical observables [26,29]. 
In particular, one can obtain an expression for the ^-dependence of the vacuum 
energy E{9). At leading order in 9, 



E{6) = le 2 X (0) = -^9 2 Jim / d'xe^ (^GG(x), ^GG(O)) , (2.21) 

where x(0) is known as the topological susceptibility. Making use of the 
anomaly relation and assuming ^> m n , this may be evaluated as [29] 



E {9) = ~9 2 m,(qq) + 0(9\ml), (2.22) 

where (qq) is the quark vacuum condensate and is the reduced quark mass, 
given by 



m u + m d 

in two-flavour QCD. This dependence on the reduced quark mass can be 
straightforwardly understood on recalling that 9 becomes unphysical as soon 
as any quark eigenstate becomes massless. Indeed, we see that the result is 
essentially fixed up to an order-one coefficient by the dictates of the anomalous 
Ward identity. On general grounds, we would expect E{9) ~ 6> 2 m*A^ ad , where 
A had is the characteristic hadronic scale required on dimensional grounds, 
which the calculation above identifies with the quark condensate. 

The quadratic dependence of the vacuum energy on 9, since its determined 
effectively by a two-point function, implies that generic CP-odd observables 
will inherit a leading linear dependence on 9. In particular, although we will 
discuss a more detailed calculation in the next section, we can obtain a similar 
order of magnitude estimate for the neutron EDM, 



d n ~ e ^ ~ 9 ■ (6 x 1(T 17 ) e cm, (2.24) 

^had 

where we identified A had = m n and used conventional values for the light quark 
masses. The experimental bound then translates into the limit, 



|0| < 10" 9 . (2.25) 

This remarkable degree of tuning in the value of 9 then constitutes the strong 
CP problem. It is aggravated by the fact that 9 is a dimensionless parameter, 
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and thus can receive corrections from unspecified sources of CP violation at 
an arbitrarily high scale. 

If we discard the possibility that this tuning is simply accidental, and search 
for a theoretical explanation for why 9 is very small, then we find that the 
existing theoretical attempts to solve the strong CP problem can be divided 
into those that are based either on continuous symmetries or on spontaneously 
broken discrete symmetries. To some extent, these two possibilities can also 
be motivated by two extreme reference points, namely when 9 is either fully 
rotated to sit in front of GG, or to manifest itself as an overall phase of the 
quark mass matrix. Although inherently basis-dependent, the former view- 
point suggests that 9 is essentially tied to the gluonic structure of QCD, while 
the latter emphasizes instead its links to the flavour sector. 



2. 3. 1 Dynamical relaxation of 9 

The energy of the QCD vacuum as a function of 9 (2.22) has a minimum 
at 9 — 0. Thus the relaxation of the ^-parameter to zero is possible if one 
promotes it to a dynamical field, called the axion [31-33]. This is motivated 
by the assumption that the Standard Model, augmented by appropriate ad- 
ditional fields, admits a chiral symmetry U(1)pq, acting on states charged 
under SU(3) C . When this symmetry is spontaneously broken at a necessarily 
high scale f a , a pseudoscalar Goldstone boson - the axion - survives as the 
only low energy manifestation. Symmetry dictates that the essential compo- 
nents of the axion Lagrangian are very simple, 



C a = l -d,ad»a + ^^GG, (2.26) 

2 Ja OH 

leading to a field-dependent shift of 9, 



9^9 + y. (2.27) 

J a 

If the effects of non-perturbative QCD are ignored, this Lagrangian possesses 
a symmetry, a — > a+const, and a is a massless field with derivative couplings 
to the SM fields, i.e. d^aip^^ip, that are not important for the solution of 
the strong CP problem. 

Below the QCD scale, one finds that U(1)pq is explicitly broken by the chi- 
ral anomaly, and thus the axion is in reality a pseudo-Goldstone boson and 
acquires a potential. The form of this potential can be read directly from 
our earlier discussion of the ^-dependence of the vacuum energy, namely 
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E{0) ~ x(0)^ 2 / 2 + •••• Accounting for the shift in (2.27), the effective ax- 
ion Lagrangian becomes, 

Cf = - 1 X (0) (W £) ■ + •••. (2-28) 

We see from (2.27) that the vacuum expectation value of the axion field (a) 
renormalizes the value of so that all observables depend on the (0 + (a) / f a ) 
combination. At the same time, such a combination must vanish in the vac- 
uum as it minimizes the value of the axion potential in (2.28). This dynamical 
relaxation then solves the strong CP problem. This cancellation mechanism 
works independently of the "initial" value of 0, which is why it is very ap- 
pealing. However, the excitations around (a) correspond to a massive axion 
particle with 



m a ~-\ X m l/ \ (2-29) 

J a 

a formula analogous to that discussed earlier for rj' (2.20). For large f a the 
axion is very light and thus has significant phenomenological consequences. 
Indeed the negative results of direct and indirect searches for "invisible axions" 
[34,35], where f a is a free scale as we have been discussing here, has now 
imposed a rather large bound, f a > 10 10 GeV, while cosmological constraints 
imply, on the contrary, that it cannot be too much larger than this (see e.g. 
[36]). 

An aspect of the axion mechanism that is perhaps not stressed as often as it 
should be is that there can be other contributions to the axion potential which 
shift its minimum away from (0 + a/f a ) = [37]. In particular, in the analysis 
above, we included only the leading term corresponding to the vacuum energy. 
However, if there are other CP-odd operators Ocp present at low scales, QCD 
effects may also generate terms linear in via nonzero mixed correlators of 
the form 



Xo 



CP 



(0) = -i lim / d 4 xe lk - x (0\T(GG(x), O CP (0))|0). (2.30) 



An example of this type is the quark chromoelectric dipole moment, Oqp = 
d q qGa'-y 5 q, appearing in (2.16). The axion potential is then modified, 



cf = - xo CP (o) (e + j) - lx(o) (e + j) + • • • , (2.31) 
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and exhibits a minimum shifted from zero. The size of this induced contribution 
to 6, i.e. #i n d = — Xe>cp(0)/x(0)> is linearly related to the coefficient of the CP- 
odd operator Oqp generating Xc C p(0)- These effects therefore need to be taken 
into account in computing the observable consequences of CP-odd sources in 
axion scenarios, and will be important for us later on. 

Before moving on, it is worth recalling that, were it realized, the simplest so- 
lution to the strong CP problem would fall into the class we are discussing, 
namely the possibility that m u = in the Standard Model Lagrangian normal- 
ized at a high scale M, or more generically, det(F u (M)) = 0. In this situation, 
the Lagrangian already possesses the appropriate chiral symmetry without the 
addition of extra fields and, as we have discussed, 6(M) then becomes unphys- 
ical. Since the identification of light quark masses is indirect, using meson and 
baryon spectra and chiral perturbation theory, the possibility that m u = has 
been debated at length in the literature [38,39], but is strongly disfavoured by 
conventional chiral perturbation theory analysis, with recent results implying 
mu/md — 0.553 ± 0.043 [40], and this conclusion is beginning to be backed up 
by unquenched (but chirally extrapolated) lattice simulations which suggest 
similar values, m u /m d = 0.43 ± 0.1 [41]. 

2.3.2 Engineering 9 ~ 

Another way to approach the strong CP problem is to assume that either P or 
CP or both are exact symmetries of nature at some high-energy scale. Then 
one can declare that 9GG be zero at this high scale as a result of symmetry. Of 
course, to account for the parity- and CP- violation observed in the SM, one 
has to assume that these symmetries are spontaneously broken at a particular 
scale Ap(cp). 

The model building problem that this sets up - one which has been made 
particularly manifest by the consistency of the recently observed CP- violation 
in P-meson decays with the KM mechanism - is that one needs to ensure that 
the subsequent corrections to 9 are small, while still allowing for an order one 
KM phase. Symmetry breaking at A P ( C p) ma Y generate the ^-term at tree 
level through e.g. imaginary corrections to the quark mass matrices M u and 
M d , 

9 ~ Arg Det(M u M d ) + ■■■ (2.32) 

where such corrections could affect either the Yukawa couplings Y u ^) or the 
Higgs vacuum expectation values, v u ( d ) (i n the SM v u — v d ), while the ellipsis 
denotes the phases of other coloured fermions. For comparison, the SM CKM- 
type phase (in basis-invariant form) is [42] 
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#km ~ Arg Det \M u Ml M d M ( 



(2.33) 



and one is then led to consider models for flavour in which the second phase 
(2.33) can be large, as is required, while the first (2.32) vanishes, or is at least 
highly suppressed. 



One class of models uses exact parity symmetry at some high-energy scale 
which implies L <-> R reflection symmetry in the Yukawa sector and thus 
hermitian Yukawa matrices which do not contribute to 9 [43]. However, this 
necessitates the extension of the SM gauge group to incorporate SU(2)r, and 
the reality of v u ^) comes as an additional constraint on the model which can 
be achieved e.g. in its super symmetric versions [44,45]. One can instead just 
demand that CP be an exact symmetry at high scales, which is then bro- 
ken spontaneously and CP- violation enters via complex vacuum expectation 
values of additional scalar fields. Models of this type can be constructed in 
which the mass matrices are complex, but have a real determinant [46,47], 
although often it is difficult to obtain a sufficiently large CKM phase. An in- 
teresting recent suggestion for getting round this problem is to use low-scale 
supersymmetry breaking [48] (see also earlier ideas [49,50]), while CP is bro- 
ken spontaneously at a much higher scale where SUSY is still exact. Strong 
interactions in the CP-breaking sector can then generate a large CKM phase, 
while a SUSY nonrenormalization theorem ensures that 6 is not generated 
until the much lower scale where SUSY is broken. 



All models of the type discussed above, that attempt to solve the strong CP 
problem by postulating exact parity or CP at high scales, have to cope with 
the very tight bound on 0. Indeed, it is not enough to obtain 6 = at tree level, 
as loop effects at and below Ap^'P) can lead to a substantial renormalization 
of the 6*-term (see, e.g. [51,52]). If the effective theory reduces to the SM below 
the scale A P ( C p), the residual low-scale corrections to the ^-term can only come 
via the Kobayashi-Maskawa phase and the resulting value for 0(5 KM ) is small. 
However, this does not guarantee that the threshold corrections at Ap^cp) are 
also small, as they will depend on different sources of CP- violation and do not 
have to decouple in the limit of large Apujp)- Such corrections are necessarily 
model-dependent. However, if the underlying theory is supersymmetric at the 
scale Ap(cp) and the breaking of supersymmetry occurs at a lower scale A SUSYl 
one expects the corrections to 6 to be suppressed by power(s) of the small ratio 
Asusy/A P ( C p) [48]. 

To summarize this section, we comment that the way the strong CP problem 
is resolved affects the issue of how large additional non-CKM CP violating 
sources can be. The axion solution, as well as m u = 0, generically allows for 
the presence of arbitrarily large CP violating sources above a certain energy 
scale. This scale is determined by comparison of higher-dimension CP-odd 
operators (i.e. dim> 5) induced by these sources with the current EDM con- 
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straints. On the contrary, models using a discrete symmetry solution to the 
strong CP problem usually have tight restrictions on the amount of additional 
CP- violation even at higher scales in order to avoid potentially dangerous con- 
tributions to the 0-term. 



3 QCD Calculation of EDMs 

Having discussed the #-term in detail, we now take a more general approach 
and consider all the relevant operators up to dimension six in the CP-odd 
Lagrangian (2.13), and move to the next level in energy scale in Fig. 1. To 
proceed, we need to determine the dependence of the nucleon EDMs, pion- 
nucleon couplings, etc., on these quark-gluon Wilson coefficients normalized 
at 1 GeV, i.e. 

d n = d n (9, di, di, w, Cij), 

9wNN = 9wNn(0, d q , w, Cij), (3.34) 

The systematic project of deducing this dependence was first initiated some 
20 years ago by Khriplovich and his collaborators, and is clearly a nontrivial 
task as it involves nonperturbative QCD physics. It is nonetheless crucial in 
terms of extracting constraints, and in particular one would like to do much 
better than order of magnitude estimates so that the different dependencies 
of the observable EDMs may best be utilized in constraining models for new 
physics. 

It is this problem that we will turn to next. In order to be concrete, we will limit 
our discussion to the nucleon EDMs and pion-nucleon couplings. The electron 
nucleon couplings, of which C s plays the most important role for the EDM of 
paramagnetic atoms, receive contributions from the semi-leptonic four-fermion 
couplings C qe in (2.18), which may be determined straightforwardly using low- 
energy theorems for the matrix elements of quark bilinears in the nucleon (See 
e.g. Refs. [53]). 

Before we delve into some of the details of these calculations, it is worth 
outlining a checklist of attributes against which we can compare the various 
techniques available for these calculations. We list below several features that 
such techniques would ideally possess: 

• Chiral invariance, including the relevant anomalous Ward identities, pro- 
vides a very strong contraint on the manner in which CP-odd sources may 
lead to physical observables in the QCD sector. As an example, distributng 6 
arbitrarily between GG and qij^q cannot alter the prediction for d n {6), and 



19 



the answer must also depend on the correct combination of quark masses, 
namely the reduced mass m*. These symmetry constraints are therefore very 
powerful, and allow a consistency check of the QCD estimates. Calculations 
are therefore more transparent if these constraints can be "built in" . 
In addition to these chrial properties, the need to deal first of all with the 
tuning of means that ideally the procedure should also corrrectly account 
for additional contributions generated under PQ relaxation. As argued in the 
previous section, the presence of CP-odd operators can shift the position of 
the axion expectation value, leading to a new class of contributions, d n (9 ind ). 
More specifically, this is the case for the CEDM sources, the presence of 
which implies that 6 must be substituted not by zero but by 6> in d given by 
[54]: 

Sind—f E ^ (3-35) 

_ ill in 

q=u,d,s y 

independently of the specific details of the axion mechanism. Here mg de- 
termines the strength of the following mixed quark-gluon condensate, 

g s (qGaq) =ml{qq), (3.36) 

Thus PQ symmetry may lead to additional vacuum contributions to the 
EDM. 

In order to make consistent use of the EDM constraints on fundamental 
CP-odd phases, it would be desirable to have the same method available 
for obtaining estimates of d n and g n NN i n terms of the relevant Wilson 
coefficients. Since different techniques have different sources of errors, use 
of the same method may allow a reduction in the uncertainty between the 
relative coefficients which, given a suite of different constraints, is ultimately 
more important then the overall uncertainty. 

Ideally, the method used should allow for a systematic estimate of the preci- 
sion of the relevant QCD matrix elements, i.e. within a framework allowing 
for a treatment of higher order subleading corrections. 
The generic dependence of d q on m q poses an additional challenge for 
obtaining precise results, through the poorly known values of the light 
quark masses, and their strong dependence on the QCD normalization 
scale. This uncertainty can be ameliorated given a method which gener- 
ates answers which depend only on scale invariant combinations such as i. e. 
(m u + m d )(qq) = -flm\ or m u /m d . 

For EDMs, the contribution of operators that are not suppressed by the light 
quark masses, m u and m^, is of considerable phenomenological interest. An 
ideal method would lead to a quantitative prediction for whether d s or w can 
compete with the contributions of light quark EDMs and CEDMs. A well- 
known example in the + channel suggests that this may indeed happen: 
i.e. {N\m s ss\N) > {N\m q qq\N) , where q = u,d. 
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Since this is a nonperturbative QCD problem, the tools at our disposal are lim- 
ited. Ultimately, the lattice may provide the most systematic treatment, but 
for the moment we are limited to various appoximate methods and none that 
are currently available can satisfy all of the demands listed above, although 
we will argue that combining QCD sum-rules with chiral techniques can sat- 
isfy most of them. While one can make use of various models of the infrared 
regime of QCD, we prefer here to limit our discussion to three (essentially) 
model-independent approaches, which vary both in their level of QCD input, 
and in genericity as regards the calculations to which they may be applied. 

However, we will first recall what is perhaps the most widely used approach for 
estimating the contribution of quark EDMs to the EDM of the neutron. This 
is the use of the SU(6) quark model, wherein one associates a nonrelativistic 
wavefunction to the neutron which includes three constituent quarks and al- 
lows for the two spin states of each. Obtaining the contribution of quark EDMs 
to d n then amounts to evaluating the relevant Clebsch-Gordan coefficients and 
one finds, 

d n {d q )^ M = l -{Ad d -d u ). (3.37) 

Although one may raise many questions regarding the reliability, and expected 
precision, of this result, we will emphasize here only the significant disadvan- 
tage that this approach cannot be used for a wider class of CP-odd sources, 
relevant to the generation of d n and g n NN- 

3.1 Naive dimensional analysis 

Although historically not the first, conceptually the simplest approach is a 
form of QCD power-counting which goes under the rather unassuming name 
of "naive dimensional analysis" (NDA) [55]. This is a scheme for estimat- 
ing the size of some induced operator by matching loop corrections to the 
tree level term at the specific scale where the interactions become strong. In 
practice, one uses a dimensionful scale A had ~ 471"/^ characteristic of chiral 
symmetry breaking, and a dimensionless coupling Ahad/A to parametrize the 
coefficients. The claim is that, to within an order of magnitude, the dimen- 
sionless "reduced coupling" of an operator below the scale A had is given by the 
product of the reduced couplings of the operators in the effective Lagrangian 
above A had which are required to generate it. The reduced couplings are de- 
termined by demanding that loop corrections match the tree level terms, and 
for the coefficient cq of an operator O of dimension D, containing N fields, 
is given by (47r) 2_JV A^J 4 C0. A crucial, and often rather delicate, point is the 
precise scale at which one should perform this matching. Within the quark 
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sector, the identification of this scale with Anf^ often seems to work quite 
well. However, for gluonic operators, the implied matching occurs at a very 
low scale where g s is very large, up to g s ~ An, and NDA has proved more 
problematic in this sector. 

To illustrate this approach, let us consider the neutron EDM induced by 9, 
in this case realized as an overall phase 9 q of the quark mass matrix, and 
also the EDM and CEDM of a light quark. The dimension five neutron EDM 
operator has reduced coupling <i n Ahad/(47r). Above the scale Ahad we need the 
reduced couplings of the electromagnetic coupling of the quark, e/(4-7r), and 
the CP-odd quark mass term, 6 q m q / A had . Thus we find, 



d n (^,/i)~e^(/i)^^, (3.38) 

A had 

where the //-dependence reflects the choice of matching scale. To obtain a 
similar estimate for the contribution of a light quark EDM, we note simply 
that it has a reduced coupling given by c^Ahad/ (47r) and thus 



dn{d q ,n) ~ d q (n), (3.39) 

which can be contrasted with the quark model estimate above. The contri- 
bution of the quark CEDM is similar, but one needs in addition the reduced 
electromagnetic coupling of the quark, e/(47r), so that 



^K,/i) = ^^<(/i), (3.40) 



where we have redefined the CEDM operator so that d q = g s d' q . This makes 
the factor g s explicit, which seems crucial to the success of NDA for gluonic 
operators as the matching needs to be performed at a large value of g s , e.g. 
g s ~ 47r as noted above. 

These examples indicate on one hand the simplicity of this approach and also 
its general applicability, but also the fact that it does not easily allow one 
to combine different contributions into a single result for the neutron EDM. 
In particular, these estimates have uncertain signs and thus can only be used 
independently with an assumption that the physics which generates them does 
not introduce any correlations. This will not generically be the case. 
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Fig. 2. Chirally enhanced contribution to the neutron EDM.. 
3.2 Chiral techniques 

Historically, the first model-independent calculation of the neutron EDM [56] 
made use of chiral techniques to isolate an infrared log-divergent contribution 
in the chiral limit (for an earlier bag model estimate see [57]). This was one 
of the landmark calculations which made the strong CP problem, and indeed 
the magnitude of the required tuning of 9, quite manifest. 

The basic observation was that, given a CP-odd pion-nucleon coupling g n NN, 
one could generate a contribution to the neutron EDM via a 7r~-loop (see 
Fig. 2) which was infarared divergent in the chiral limit. In reality this log- 
divergence is cut off by the finite pion mass, and one obtains, 

d ^ = ^Mj^^ <3A1) 

where A is the relevant UV cutoff, i.e. A = m p or M n . One can argue that 
such a contribution cannot be systematically cancelled by other, infrared finite, 
pieces and thus the bound one obtains on g^NN in this way is reliable in real- 
world QCD. 

This reduces the problem to one of computing the relevant CP-odd pion- 
nucleon couplings. For a given CP-odd source Oqp, we have 

(Nir a \O CP \N') = j-(N\[0 OP , Jg 5 ]\N') + rescattering, (3.42) 

JIT 

justified by the small t-channel pion momentum. The possible rescattering 
corrections will be discussed below. If we now specialize to the #-term, as 
in [56], with Oqp = —9 q m* J2f QfijsQf then the commutator reduces to the 
triplet nucleon sigma term, and we find 

= ^(Plrtb) fl " 4V (3.43) 
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Fig. 3. Two classes of diagrams contributing to the CP-odd pion-nucleon coupling 
constant. 

One can then determine (N\qr a q\N) from lattice calculations or, as was done 
in [56], by using global SU(3) symmetry to relate it to measured splittings in 
the baryon octet. 

The final factor on the right hand side of (3.43) reflects the vanishing of the 
result in the limit that the chiral anomaly switches off and 77 (or rf in the three- 
flavour case) is a genuine Goldtone mode. This factor is numerically close to 
one and was ignored in [56]. It arises because in (3.42) we should also take into 
account the fact that the CP-odd mass term can produce rj from the vacuum 
and thus, in addition to the PCAC commutator, there are rescattering graphs 
with 77 produced from the vacuum and then coupling to the nucleon, and the 
soft pion radiated via the CP-even pion-nucleon coupling [58]. 

Although this technique is not universally applicable, one can also contemplate 
computing the contribution of certain other sources, e.g. the quark CEDMs. 
Using the same PCAC-type reduction of the pion in (Nir a \Ocp\N') as in 
(3.42), one can reduce the calculation of g^NN to the matrix elements of di- 
mension five CP-even operators. In doing this, one has to take into account a 
subtlety for CEDM sources, first pointed out in [59,58,60], namely that a sec- 
ond class of contributions, the pion-pole diagrams (Fig. 3b), now contribute at 
the same order in chiral perturbation theory. In an alternative but physically 
equivalent approach, one can perform a chiral rotation in the Lagrangian to set 
(0\Cgp\n°) = 0, thus making this additional source of CP-violation explicit 
at the level of the Lagrangian [61]. 

After this consistent PCAC reduction of the pion, the intermediate result for 
9-kNN takes the following form, 

— — (N\d u (ug s Gau — rn^uu) — dd(dg s (Ga)d — rn^dd)\N) 



The second line in this expression contains the same matrix element as (3.43), 
and vanishes when PQ symmetry sets the axion minimum to 9 in & (3.35). The 
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+^(e-6 ind )(N\qr 3 q\N). 



(3.44) 
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remaining terms are proportional to the following combination, 



(N\qg s Gaq-m 2 qq\N), 



(3.45) 



which unfortunately cannot be estimated using chiral techniques, and requires 
genuine QCD input. A naive vacuum saturation hypothesis in (3.45) leads to 
the vanishing of this expression. This is a rather fundamental problem which 
limits the precision of various approaches, e.g. those based on the use of low 
energy theorems to estimate (3.45) [62,63], to obtain the dependence of g n NN 
on the CEDMs. 

This limited applicability is one problem that currently afflicts the chiral ap- 
proach. A more profound issue is that the terms enhanced by the chiral log, 
while conceptually distinct, are not necessarily numerically dominant. Indeed 
there are infrared finite corrections to (3.41) which, while clearly subleading 
for — > 0, are not obviously so in the physical regime. This dependence on 
threshold corrections has been observed to provide a considerable source of 
uncertainty [64] (see also [65]). 

3.3 QCD sum-rules techniques 

An alternative to considering the chiral regime directly, is to first start at 
high energies, making use of the operator product expansion, and attempt 
to construct QCD sum rules [66] for the nucleon EDMs, or the CP-odd pion 
nucleon couplings. This approach in principle allows for a systematic treatment 
of all the sources, and is motivated in part by the success of such approaches to 
the calculation of baryon masses [67] and magnetic moments [68] . For a recent 
review of some aspects of the application of QCD sum rules to nucleons, see 
e.g. Ref. [69]. 

The basic idea is familiar from other sum-rules applications. One considers the 
two-point correlator of currents, t)n(x), with quantum numbers of the nucleon 
in question in a background with nonzero CP-odd sources, an electromagnetic 
field F^, and also a soft pion field 7r a , 



where Q 2 = —p 2 , with p the current momentum. It is implicit here that the 
soft pion field admits PCAC reduction, and then in the case of CEDM sources 
corresponds to an external field coupled to the operator qg s Gaq — m^qq, as in 




(3.46) 



(3.44-3.45). 
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One then computes the correlator at large Q 2 using the operator product 
expansion (OPE), generalised to incorporate condensates of the fields, and 
then matches this to a phenomenological parametrization corresponding to an 
expansion of the nucleon propagator to linear order in the background field 
and CP-odd sources, and corresponding higher excited states in the relevant 
channel. In practice, one makes use of a Borel transform to suppress the con- 
tribution of excited states, and then checks for a stability domain in Q 2 , or 
rather the corresponding Borel mass M, where the two asymptotics may be 
matched. 

Let us discuss this procedure in a little more detail. For the neutron, there are 
two currents with lowest dimension that are commonly used in QCD sum-rules 
and lattice calculations, 

T)i = 2e abc (d^C-f 5 u b )d c , 

V2 = 2e abc (d^Cu b )-/ b d c , (3.47) 

of which the second, 7)2, vanishes in the nonrelativistic limit, and lattice simu- 
lations have shown that 7)1 indeed provides the dominant projection onto the 
neutron state [69,70]. In the truncated OPE expansion, an admixture of 7)2 
can nonetheless be used to optimise convergence, and thus it is natural to 
parametrize the current in the form r\ n = 7)1 + 07}2 introducing an unphysical 
parameter f3. The truncated OPE will then inherit a dependence on this pa- 
rameter, which can then be fixed to improve convergence once the sum-rule 
has been constructed for a given physical quantity. 

The correlator (3.46) exhibits various Lorentz structures (LS) in its OPE and, 
in selecting one to consider, one needs to be aware that in a CP- violating 
background the coupling of the current to the neutron state, described by a 
spinor v, is not invariant under chiral rotations, i.e. (Q\r) n |0) = \e ia ^/ 2 v. It 
turns out that of the terms contributing to the neutron EDM, there is a unique 
structure which is invariant under chiral rotations, namely LS={Fo"7 5 ,^}, thus 
this is the natural quantity on which to focus in constructing a sum rule for 
the EDM (for alternatives, see [71,72]). Correspondingly, LS=^ is the relevant 
chiral- invariant structure for the g n NN sum rule. 

In constructing a phenomenological model for the current correlator, it is 
apparent that in expanding to linear order in the external field we are effec- 
tively considering a three-point function. It is then not particularly useful 
to work with the spectral density, as is standard for two-point functions, 
and the conventional approach is to parametrize the correlator itself, i.e. 
-Qphen _ Lgy(p2^ _|_ . . . _ xhe function f(p 2 ) will in general have an expansion 
in double and single pole terms, and then a continuum modelling the tran- 
sitions between excited states. A Borel transform can be applied to suppress 
the contribution of excited states. However, a well-known complication [68] of 
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Vn(0) 



Fig. 4. A leading contribution to the neutron EDM within QCD sum rules. Sensi- 
tivity to the CP- violating source enters through the two soft quark lines which lead 
to a dependence on the chiral condensate. 

baryon sum rules in external fields is that the single pole terms, corresponding 
to transitions between the neutron and excited states, are not exponentially 
suppressed by the Borel transform and thus provide the leading contribution 
from the excited states, with a coefficient which is not sign definite. This must 
then be treated as a phenomenological parameter to be determined from the 
sum rules themselves. In this approximation, we then find [73,58,74], 



where the constants A, A' parametrise the single-pole contributions. One can 
then go further and construct a full continuum model to match the high-Q 2 
asymptotics, but as discussed below this refinement has minimal impact in 
comparison to the single pole terms A and A'. We now turn to the calculation 
of the OPE for d n and g n NN- 

• Nucleon EDM calculations 



The OPE for d n is conveniently constructed in practice by first computing the 
generalized quark propagator, expanded in the presence of the background 
field, the CP-odd sources, and also the vacuum condensates. One then com- 
putes the relevant contractions in (3.46) to obtain the OPE to the appropriate 
order. Although it would take us too far afield to describe this procedure in 
detail, we can exhibit some of the dominant physics by looking at just one 
class of diagrams which arise in evaluating the OPE for (3.46). In particular, 




(3.49) 



(3.48) 
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in Fig. 4 two of the quarks in the nucleon current propagate without interfer- 
ence, carrying the large current momentum, while the third is taken to be soft 
and so induces a dependence on the appropriate chiral quark condensates. We 
may then make use of standard arguments [29] used to determine E(9) (2.22), 
which utilize the anomaly in the axial current and the fact that m, n ^> m n , to 
determine the dependence of these condensates on the CP-odd source. 

Let us consider the contribution of 9. For additional control over the chiral 
transformation properties of the answer, we split 9 into several terms in the 
Lagrangian, 

^dim=4 = G G^G^ a - ]T mM^i (3.50) 

i=u,d,s 

The diagram in Fig. 4 then leads to a dependence on 



™<?(O|gcV75g|O)0,F = im*9 G (0\qa^q\0) F + 0{ml) 

= ix q G m*F^{qq) + 0{ml). (3.51) 

In the first equality the dependence on 9q has been determined as in the 
standard computation of E(9) in (2.22) [29], where the terms of 0(ml) are 
subleading only because there is no U(l)-axial Goldstone mode, i.e. m T] m n . 
The dependence of (3.51) on 9q rather than 9 means that there are additional 
contributions from 9 q that have to be taken into account, at this order in the 
OPE, in order to restore the dependence on the physical 9 combination. In 
fact, it turns out that the procedure is generally more complicated and one 
must also incorporate the the mixing of the currents rji and i]2 (3.47) with their 
chirally-rotated (or CP-conjugate) counterparts [73,58] in order to restore the 
dependence on 9 and exclude unphysical parameters such as 9 G — 9 q . 

In the second line of (3.51) we have introduced the so-called electromagnetic 
susceptibility of the quark vacuum, Xq, given by 



(qcr^F = XgF^(qq), (3.52) 

and chosen for simplicity to be favour independent, Xq = X e q [68]. In fact, Xq 
is one among a whole series of mixed quark-gluon condensates that have to be 
taken into account in the calculation of the OPE for d n . However, it turns out 
that numerically x is rather large, x — ~ (5 — 9) GeV~ 2 [75,76], and results 
in the diagram of Fig. 4 being numerically very important. The remaining 
condensates are numerically quite small in comparison, and we will neglect 
them in what follows, referring the reader to Refs. [77,73,58,60,74] for more of 
the details involved in these calculations. 
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The leading order and next-to-leading order contributions to the OPE induced 
by the CP-odd operators of dimension four and five are given by 



% 



n OPE (g 2 ) = -^(gg){^75,^} 



x ( ln( -v 2 ) + - j 2 In (-^ ) [tt nlo ] ) . (3.53) 



At leading order the quark EDMs induce tt^q 



10 + 6ft 



- d h 



3 + 2(3 - ft 



(3.54) 



while the 0-term, and the CEDMs, are responsible for ttj*q, 



Trgj =4(1 + ft\ d m d P d - (1 + ft\ u m u P u 

+2(1 - ft)m^xu + Xd)(P u - P d ), (3.55) 

where P q = 9 q + ^^ c,p contains the dependence on the CP-odd sources. 

The next-to-leading order terms in (3.53), 7Tnlo; are associated with a class 
of diagrams in which one of the propagators exhibits a logarithmic infrared 
divergence as can be seen in (3.53). The magnitude of these terms is clearly 
ambiguous through the logarithmic dependence on the cutoff A IR , although in 
practice with a cutoff A IR ~ Aq C d the logarithms are not particularly large 
in the momentum regime for which the resulting sum rules are optimal. We 
will not need to exhibit them explicitly here, but it is important that for all 
of these terms 7Tnlo oc (1 — (3), which differs from the (1 + ft 2 dependence of 
many of the dominant leading order terms. 

The truncated OPE in (3.53) necessarily depends on the unphysical parameter 
(3, which can therefore be chosen in such a way as to optimise the convergence 
of the expansion. Such optimization problems arise in many areas of physics 
[78], and in the present case with only two terms of the series in hand, the 
most practical approach is to use "fastest apparent convergence" (FAC), which 
involves choosing (3 to set the highest known term in the series to zero. Histor- 
ically, Ioffe [79] introduced an FAC-like criterion in analysis of the mass-sum 
rule, which has proved quite successful. For the CP-even sector this involves 
the choice (3 — — 1 to cancel the subleading terms. We shall also follow this 
approach as it has the added advantage of canceling the ambiguous infrared 
logarithmic terms. As first discussed in [73,58], for the CP-odd sector this in- 
volves the choice (3 = 1. The difference in (3, as compared with the choice in the 
CP-even sector, is not surprising as (3 is unphysical and there is then no reason 
to expect optimal choices to be the same for different physical observables. 
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A remarkable feature of the choice (3 = 1 becomes apparent if we now rewrite 
the OPE expression for (3 = 1. All the subleading logarithmic terms in (3.53) 
are now cancelled, while the leading order terms adopt the elegant form, 



nftli = 4 [4xdm d P d - Xum u P u ] 

nl =1 =A[Ad d -d u ]. (3.56) 

It is remarkable that the contribution of the u quark in each term is precisely 
— 1/4 that of the d quark, which is the combination suggested by the SU(6) 
quark model! One might suspect that this is due to the minimal valence quark 
content of the current, but the fact that this structure arises only for the 
choice (3 = 1 is rather surprising given that only the 771 current survives in the 
nonrelativistic limit. 

One of the most important aspects of the whole calculation is the consis- 
tent treatment of the CP-odd vacuum condensates entering via the CP-odd 
sources P q . This calculation can be done relatively easily with the use of chi- 
ral techniques, and at leading order in mh^/m^,. A useful constraint on this 
calculation is that the anomalous chiral Ward identity must be respected for 
each quark flavour. This provides a useful check if we consider artificially de- 
coupling the s-quark, while at the same time sending either m u or m d to zero. 
In this regime, the dependence on 9 in particular must vanish, which fixes the 
remaining quark mass dependence in terms of m*. The final result for the P u 
and P d sources, in the limit of m u ^ <C m s , reads 



m u{d) P u(d) - m*V — — , (3.57) 

^ V m d(u) m s J 

which respects the anomalous chiral Ward identity. If 9 is removed by PQ 
symmetry, i.e. substituted by 9 in d (3.35), Eq. (3.57) simplifies even further, 



m pPQ _ _mMi 

" L q r q — 9 • 



(3.58) 



Putting all the ingredients together, after a Borel transform of (3.53) and 
(3.48), and using (3 = 1 as discussed earlier, we obtain the sum rule, 



X m n d n + AM' 



M 4 



32tt 2 



3 m2/M 2 /- 



(qq) 



7T 



(X) 
0=1 



+ 7T 



0=1 



+ 0(M°), (3.59) 



where M is the Borel mass, and the relevant contributions are given in (3.56). 
Clearly, the presence of three parameters d n , A and A in (3.59) necessitates 
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the use of additional sum rules, and the coupling A is conveniently obtained 
from the well-known sum rules for the two-point correlation function of the 
nucleon currents in the CP even sector (see e.g. [69] for a review). 

Rather than reviewing the full analysis, let us consider a simple estimate ob- 
tained from the leading order terms in the OPE of (3.59) a la Ioffe's derivation 
of the nucleon mass formula [67]. We set A = 0, and taking M = m n ~ 1 GeV, 
we divide the sum rule (3.59) by the standard CP-even sum rule for A ob- 
tained for the Lorentz structure j> and (3 = 1 (the choice (3 = — 1 in the latter 
sum rule leads to a similar result). The resulting estimate takes the following 
form, 



d 



'CSt 



8^\(qq}\ 



— #ind) 



-3(4^ - d u ) + —^—(Ae d d d - e u d u ) 



(3.60) 



where 6 m & again is a linear combination of d q /m q (3.35). The coefficient in 
front of the square brackets in (3.60) is very close to 1, given Ioffe's estimate 
for m n , rr? n ~ 87r 2 |(gg)| [67]. Indeed, this estimate shows no deviation at all 
from the naive quark model result for d n (d q )\ Using Vainshtein's value for 
X = —N c /(4:iT 2 f%) ~ —9 GeV -2 [76], obtained using pion-dominance for the 
longitudinal part of certain anomalous triangle diagrams, along with the Ioffc 
formula for m n , the estimate for d n {9) becomes 



d 



est 



em*9 



(3.61) 



which coincides with the chiral estimate (3.41) if 5u(£>|<? r3 Q , |p) ln(A/m n ) is of 
order 2, where qa — 9-kNN fw / m n- Needless to say that within the accuracy of 
both methods the two estimates coincide. If 6 is removed by PQ symmetry, 
then within the same approximation the resulting estimate reads 



df = ~d d 



\d u 



2mle 



m n (m u + m d ) 



(fdd + \du) , 



where the approximate relation itLq 



(3.62) 



-mz 



has been used, along with (m u + 



m d)\(qq)\ = f1 m \-> are used. One immediately sees that the CEDM contribu- 
tions are significant and comparable in magnitude in fact to the effects induced 
by quark EDMs. 



We can give a more precise numerical treatment by making use of the following 
parameter values: For the quark condensate, we take a central value of 
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(0|gg|0) = -(0.225 GeV) 3 , 



(3.63) 



while for the condensate susceptibilities, we have [75] 



X ^ -5.7 ±0.6 GeV" 2 , 

m 2 ~ -0.8 GeV 2 . (3.64) 

This determination of x is based on spectral sum-rules [75] and is slightly 
lower than the value obtained in [76]. 

A more systematic treatment of the sum-rule [73,58,74] indicates that a sta- 
bility domain exists for relatively low Borel mass scales, of M ~ 0(O.8GeV). 
Convergence of the OPE is apparently not in danger as this low scale arises 
via the two step procedure used, in which the OPE is naturally formulated 
around the neutron scale, M ~1 GeV, while the chiral techniques used to ex- 
tract the dependence of the condensates on the CP- violating sources lower the 
effective scale, but also introduce additional combinatoric suppression factors 
as the dimension of the condensates increases. In order to test the stability of 
the sum rule, and obtain an estimate for the uncertainty due to the handling 
of excited states, one can generalise the expression for (3.59), by including 
a more systematic parametrisation of the continuum, and also by including 
1-loop anomalous dimensions for the currents and condensates entering the 
sum-rule. However, as alluded to earlier, these refinements have a rather min- 
imal impact, moving the stability domain by no more than 10-15%. This is 
relatively small compared to the primary sources of error, namely the satura- 
tion hypothesis for the condensates, the need to extract the single-pole term 
from the sum-rules and, perhaps most significantly, the dependence on (5. 

Extracting a numerical central value from the sum-rule, employing numerical 
estimates for the condensates (3.64), and estimating the precision through 
consideration of the sources of error listed above, we find the results first 
presented in [73,74], 



d n (0) = (1 ± 0.5)— ^^—9 x 2.5 HT 16 ecm, (3.65) 



(225MeV) 3 

,_!« 

(225MeV) 3 



dn Q (d q , d q ) = (1 ± 0.5) J ^fj \l.le(d d + 0.5d u ) + lA(d d - 0.25d u ) 



where we intentionally split the formula into two parts, d n {9) and d n (d q ,d q ) 
in the presesence of PQ symmetry. In the generic case, the two lines in (3.65) 
must be added together and 9 substituted by 9 — 9 ind . 

The result (3.65) offers several interesting consequences. Note that the overall 
factor of (qq) combines with the light quark masses from short-distance ex- 
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pressions for di and di to give a result ~ f 2 m 2 (l + 0{m u /m<i)) thus reducing 
the uncertainty due to poor knowledge of the absolute values of quark masses 
and condensates. The constribution from the strange quark CEDM is sup- 
pressed by the m*/m s ratio, and is completely removed at this order once the 
PQ symmetry is switched on. However, the use of nucleon currents with only 
valence quarks leads one to suspect that gluon and sea-quark contributions 
may be under-estimated, since they enter only at higher orders. It is possible 
that the question of d n (d s ) may be resolved in future lattice simulations, given 
an appropriate lattice implementation of chiral symmetry. 

Compared to the techniques outlined previously, this approach has the signif- 
icant advantage that all of the sources up to dimension five can be handled 
systematically and thus relative signs and magnitudes can be consistently 
tracked. As indicated in (3.65), one can also make a systematic estimate of 
the precision of the result, where the errors are due to the contribution of ex- 
cited states, neglected higher dimensional operators in the OPE, and also an 
ambiguity in the nucleon current. Comparing the numerical result with those 
obtained using NDA and chiral techniques one finds, as is to be expected, that 
the results agree in terms of order of magnitude. In particular, our result for 
d n implies a bound \9\ < 3 x 1CT 10 . 

In progressing to consider the contribution from sources of higher dimension, 
problems arise through the appearance of certain infrared divergences at low 
orders in the OPE, while a number of unknown condensates also enter and 
render a corresponding calculation for dimension six sources intractable. One 
can nonetheless estimate the contribution of these operators by utilizing a 
trick which involves relating the EDM contribution to the measured anomalous 
magnetic moment /i n via the 75-rotation of the nucleon wavefunction induced 
by the CP-odd source [37], 



(N\6C CP \N) 

dn ~ ^ m n Ni l5 N ■ (3 - 66) 

One may analyze the 75-rotation using conventional sum-rules techniques, and 
for the Weinberg operator, one can obtain the following estimate [80] 



KM I ~ \fi n \^0w\n(M 2 /fi 2 IR ) ~ e 22 MeV w(l GeV), (3.67) 

taking M//j, IR = 2, where M is the Borel mass and fj,m is an infrared cutoff, 
and g s = 2.1. 

We can also apply this technique for the contribution of four-fermion oper- 
ators. For SUSY models with generic parameters CP-odd four-fermion op- 
erators are negligible due to the double helicity-flip requiring an m 2 depen- 
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dence and rendering these operators effectively of dimension eight. However, 
for large tan (3, there are enhancements for operators proportional to with 
i,j = d,s,b which can partially overcome this suppression thus altering the 
conventional picture of EDM sources (see Fig. 1). An important class of con- 
tributions in this case involves the four-ferimon operators with a 6-quark. The 
contribution of these sources to d n can again be estimated using the same 
technique as above [81], 

K(C y )| ~ e 2.6 x l(T 3 GeV 2 + .75^^) . (3.68) 

We should emphasize that both the dimension six estimates above necessarily 
have a precision not better than (9(100%), and one cannot reliably extract the 
sign. Fortunately, the numerical size of these dimension six contributions is 
often negligible, and thus does not significantly impact the phenomenological 
application of EDM constraints. 

• Calculation of g n NN 

The other primary source of nuclear EDMs, leading to the observable EDMs 
in diamagnetic atoms, arises through nucleon interactions mediated by pion- 
exchange with CP violation in the pion- nucleon vertex. As discussed in the 
preceding subsection, the calculation of these couplings involves essentially 
two steps: the first is a PCAC-type reduction of the pion in (Nir a \Ocp\N') 
as in (3.42), and the second is an evaluation of the resulting matrix elements 
over the nucleon. It is this second part for which QCD sum-rules may usefully 
be employed, and here we sketch its application to the computation of the 
dependence of g n NN on dimension five CP-odd sources in (2.16) [61]. 

The main difficulty, as alluded to in the previous section, is the partial cancel- 
lation between the rn\qq and qg s Gaq sources in (3.44). Within the sum-rule 
analysis, one can nonetheless trace this cancellation up to the next-to-next- 
to-leading order [61]. Indeed, following the approach outlined earlier of fixing 
the unphysical parameter (3 to suppress the highest calculated order in the 
OPE, leads again to (3 — 1. However, in this case, the choice (3 = 1 sets all the 
LO, NLO, and NNLO terms to zero! This cancellation is, however, seemingly 
an artefact of the purely valence-quark content of the current, and is not dic- 
tated by symmetry. Thus, although the result will necessarily have a strong 
dependence on (3. one can obtain a numerical estimate by varying f3 through 
an appropriate range. To get some idea of the size of the result, let us work 
only with the current r\\ which survives in the nonrelativistic limit and has 
the dominant overlap with the neutron state, i.e. set (3 = 0. Assuming also 
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the dominance of the nucleon double-pole term and the leading order term in 
the OPE, and separating different isospin structures, we find [61] 



£L ~ \ j^(du - d d ), g^ N ~ | + d d ). (3.69) 

Unless d u — d~d — 0, the CP-odd coupling g^NN * s several times larger than 
gf^N- Concentrating therefore on g^NNi we have 



W lO^cm (225MeV)3 O.SGeV 2 ' [ ' 

where we re-instate the dependence on the relevant condensates, normalized 
to their central values. The estimate (3.70) is half the size of the value for 

9*nn used in t 59 ]- 

The full numerical treatment of the sum-rule for the isospin-one coupling, 
produces the following result at next-to-next-to- leading order [61], 



W - l_ x X 1U 1Q „ 26cm (225MeV)3 » ^ n ) 

where the poor precision is essentially due to the cancellation of the leading 
terms as described above. We emphasize that a more precise calculation of the 
matrix element (3.45) would significantly enhance the quality of constraints 
one could draw from the experimental bounds on diamagnetic EDMs, and 
thus constitutes a significant outstanding problem. Again, it seems progress 
may have to wait for further developments in lattice QCD. 

In considering the contribution of dimension six sources to g w NN, we note 
firstly that the three-gluon Weinberg operator GGG is additionally suppressed 
by m q and can be neglected. However, as for d n , for SUSY models with large 
tan/3, certain four-fermion operators C qiq2 may be relevant, and can be ob- 
tained via vacuum factorization, as the two diagrams in Fig. 3 now fail to 
cancel, 



g^ N ~-8x l(T 2 GeV 2 + 3.3k^ + ^(1 - 0.25k)) , (3.72) 

V m d m s m b J 

where k = (N\m s ss\N) / 220 MeV, with the preferred value k ~ 0.5 [53]. 
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Fig. 5. A particular 3-loop contribution [82] to the d-quark EDM induced by the 
KM phase in the standard model. The box vertex denotes a contacted VF-boson line 
connected to the light quarks, while it is implicit that the external photon line is to 
be attached as appropriate to any charged lines. 

4 EDMs in models of CP violation 

We have now moved to the highest level in Fig. 1, which is where the EDM 
constraints can be applied to directly constrain new sources of CP violation. In 
this section, we will breifry discuss these constraints, firstly looking at why the 
Standard Model itself provides such a small background, and then why most 
models of new physics, and supersymmetry in particular, tend to overproduce 
EDMs and are thus subject to stringent constraints. 



4-1 EDMs in the Standard Model 

The recent discovery and exploration of CP violation in the neutral P-meson 
system [7] is, along with existing data from C P- violation observed in K- 
mesons, (within current precision) in accord with the minimal model of CP 
violation known as the Kobayashi-Maskawa (KM) mechanism [3]. This intro- 
duces a 3 x 3 unitary quark mixing matrix V in the charged current sector of 
up and down-type quarks taken in the mass eigenstate basis, 



This model possesses a single CP- violating invariant in the quark sector, 
J cp = I m (VtbV t * d V c dV* b ) — 3x 1CT 5 . This combination, along with 6*qcd, are the 
only allowed sources of CP violation in the Standard Model (treating "Stan- 
dard Model neutrinos" as massless). In addition to this, CP violation in the 
SM vanishes in the limit of an equal mass for any pair of quarks of the same 
isospin, e.g. d and s, u and c, etc. These two conditions are extremely powerful 
in suppressing any KM-induced CP-odd flavour-conserving amplitude. 




(4.73) 
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Fig. 6. A leading contribution to the neutron EDM in the Standard Model, arising 
via a four-quark operator generated by a strong penguin, and then a subsequent 
enhancement via a chiral tt + loop. 

• quark and nucleon EDMs 



The necessity of four electroweak vertices requires that any diagram capable 
of inducing a quark EDM have at least two loops. Moreover, it turns out 
that all EDMs and color EDMs of quarks vanish exactly at the two-loop level 
[83], and only three-loop diagrams survive [84,82], as in Fig. 5. A leading-log 
calculation of the three-loop amplitude for the EDM of the <i-quark produces 
the following result [82], 



m d m c a s G F J C p 2( 2, 2\-i 1 2\ iaha\ 

= 6 108^ ( m b/m c ) \n(M w /m b ). (4.74) 

Upon the inclusion of the other contributions, it produces a numerical estimate 

d KM „ 1(r 34 e cm _ (475) 



The only relevant operator that is not zero at two-loop order is the Weinberg 
operator [85], but its numerical value also turns out to be extremely small. 
Indeed the largest Standard Model contributon to d n comes not from quark 
EDMs and CEDMs, but instead from a four-quark operator generated by a 
so-called "strong penguin" diagram shown in Fig. 6. This is enhanced by long 
distance effects, namely the pion loop, and it has been estimated that this 
mechanism could lead to a KM-generated EDM of the neutron of order [86], 



d KM „ 1Q -32 e cm _ (476) 

However, this is still six to seven orders of magnitude smaller than the current 
experimental limit. 
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• lepton EDMs 

The KM phase in the quark sector can induce a lepton EDM via a diagram 
with a closed quark loop, but a non-vanishing result appears first at the four- 
loop level [87] and therefore is even more suppressed, below the level of 

d KM < 10 -38 e (477) 

and so small that the EDMs of paramagnetic atoms and molecules would be 
induced more efficiently by e.g. Schiff moments and other CP-odd nuclear 
momenta. 

In this regard, we note that recent data on neutrino oscillations points toward 
the existence of neutrino masses, mixing angles, and possibly of new CKM- 
like phase(s) in the lepton sector. Under the assumption that neutrinos are 
Majorana particles, the presence of these new CP-odd phases in the lepton 
sector allows for a non- vanishing two-loop contributions to d e [88], without 
any further additions to the Standard Model. However, recent calculations 
[89] show that a typical see-saw pattern for neutrino masses and mixings only 
induces a tiny contribution to the EDMs in this way, of 0(m e mlG%), unless 
a fine-tuning of the light neutrino masses is tolerated in which case d e could 
reach 10~ 33 e cm. Therefore, within this minimal extension of the Standard 
Model allowing for massive neutrinos, the electron EDM is not the best way 
to probe CP violation in the lepton sector. 

• Probing the scale of new physics 

The Standard Model predictions for EDMs described above are well beyond 
the reach of even the most daring experimental proposals. This implies in 
turn that the Kobayashi-Maskawa phase provides a negligible background and 
thus any positive detection of an EDM would necessarily imply the presence 
of a non-KM CP- violating source. Before we consider some of the models 
which provide motivations for anticipating such a discovery, it will first be 
useful to consider in more general terms how high an energy scale one could 
indirectly probe with EDM meaurements. Indeed, we are led to ask first of 
all, what energy scale of new CP- violating physics is probed with the current 
experimental sensitivity to EDMs? Secondly, given the small KM background, 
we might also ask for the largest energy scale that could be probed in principle 
before reaching the level where the Standard Model KM contibutions would 
become significant. 

To try and answer these questions in a systematic way, let us consider a toy 
model containing a scalar field (which is Higgs-like, but needn't be the SM 
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Higgs) coupled to the SM fermions with scalar and pseudoscalar couplings, 
£ int = —J2i ( P4'i(yi + izi^ipi- The presence of both scalar and pseudoscalar 
couplings i/i and Zi breaks CP invariance. 

Assuming that the scalar mass M is large, we integrate this field out and 
match the resulting coefficients with the Wilson coefficients listed in (2.16) 
and (2.18). In particular, at tree level, we obtain the following contribution to 
(2.6), 



Cf ^ ^z e (3(y d + y u ) + k Vs + • • • ), (4.78) 

where the ellipsis stands for the contribution of heavy quark flavours, and the 
couplings are normalized at 1 GeV. If we now make the assumption that there 
is no correlation with other sources of CP violation, e.g. the electron EDM 
d e , then with the use of the experimental constraint on dn and the results 
of atomic claculations (2.8), we arrive at the following limit on the CP-odd 
combination of couplings and the mass M, 



Given the most optimistic assumption about the possible size of these cou- 
plings, i.e. z e y q ~ 0(1), we can conclude that the current experimental EDM 
sensitivity translates to a bound on M as high as M CP ~ 10 6 GeV. If instead 
we insert the largest Kobayashi-Maskawa prediction for the Tl EDM of order 
~ 10~ 35 e cm in place of the current sensitivity, we obtain Mgp x ~ 10 11 GeV, 
as the ulmitate scale which can be probed via these dimension 6 operators 
before the onset of the "KM background" . 

For comparison, allowing for arbitrarily large flavour-violating couplings of <fi 
to fermions, we can also deduce the sensitivity level to New Flavour Physics 
(NFP) in a similar way. For example, requiring that four-fermion operators 
that change flavour by two units, e.g. s^ 5 ds^ 5 d and the like, do not introduce 
new contributions to Am B , Am K and e K that are larger than the SM con- 
tribution, one typically finds that M NFP > 10 7 — 10 8 GeV in this scenario. 
Thus, we see that the sensitivity of EDMs is already approaching this bench- 
mark and, unlike the contraints from the AF = 2 sector, can be significantly 
improved in the near future. 

At this point, we should emphasize that in this example, we are relaxing all 
constraints on the flavour structure by allowing order one couplings of the 
scalar field <fi to the light fermions. These couplings violate chirality maxi- 
mally, and if they were part of a more realistic construction, for example a 
two-Higgs doublet model, one would expect that yi and z; t will have to scale 
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according to the fermion mass mj. In this case, the sensitivity to M clearly 
drops dramatically, and the tree-level interactions (4.78) are not necessarily 
the dominant contributions to EDMs, as heavy flavours may contribute in a 
more substantial way via loop effects [90,91]. Indeed, if new physics above the 
electroweak scale preserves chirality, as is often assumed, one expects that for 
the light flavours d,~ex(l - 10) MeV/M 2 . Taking the electron EDM, and 
the Tl EDM bound, as a concrete example we find under this more restrictive 
assumption that 



and consequently the current level of sensitivity to new CP- violating chirality- 
preserving physics drops somewhat, but for reference this scale is still well- 
beyond the centre-of-mass energy of the LHC. 

If we put the current EDM bounds into the broader context of precision tests 
of the Standard Model, we see that the present bounds in Table 1 imply that 
EDMs occupy an intermediate position in sensitivity to mass scales for new 
physics, between the electroweak precision tests (EWPTs) and flavour viola- 
tion in AP = 2 processes noted above. The EWPTs from LEP impose a bound, 
-^ewpt > few TeV, through constraints on various dimension six operators, 
e.g. oblique corrections to gauge boson propagators. Since Mewpt Mz, 
this has been dubbed the "little hierarchy problem" [92]. Indeed, while there 
are general expectations that the Standard Model is an effective theory, and 
will be corrected at scales of about a TeV, it is clear from this discussion that 
precision constraints in many sectors do not contain any hints of new physics 
beyond the Higgs at the weak scale, and in this sense EDMs are no excep- 
tion. The remarkably large scale Mqp implied by EDM limits requires, at least 
within our current level of understanding, a tuning in the CP-odd sector of 
physics beyond the standard model that we lack a coherent explanation for. 
The recent data from BaBar and Belle on CP violation in the neutral P-meson 
sector, which thus far is consistent with the KM model, within which CP vio- 
lation is maximal within the confines of the flavour structure, only makes this 
tuning more pronounced, since we lack a strong motivation to enforce any 
additional CP-violating phases to be small. Moreover, further experimental 
progress in the near future could, given null results, push the value of Mqp 
closer to M NFP . From this viewpoint EDMs provide our most powerful tool in 
probing the question of whether CP- violation and flavour physics are intrinsi- 
cally linked, as indeed they are within the electroweak Standard Model. This 
issue stands out as one of the most important ways in which EDMs may assist 
in demystifying some of the less constrained parts of the Standard Model. 
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4-2 EDMs in supersymmetric models 

Having demonstrated the generic importance of EDM constraints for TeV- 
scale physics, we would now like to make this analysis more concrete by fo- 
cussing on models with electroweak scale superymmetry and reviewing their 
predictions for EDMs (see e.g. [93] for reviews of EDMs within several other 
classes of models). 

Supersymmetric extensions of the Standard Model provide perhaps the most 
natural solution to the gauge hierearchy problem by automatically cancelling 
the quadratically divergent contributions to the Higgs mass. Supersymmetry 
is thought of here as a symmetry of nature at high energies, whereas at the 
electroweak scale and below it is obviously broken. Ensuring that supersym- 
metry breaking does not re-introduce quadratic divergences, and is compatible 
with the observed low energy spectrum, still allows for a large number of new 
dimensionful parameters, unfixed by any symmetry, that are usually called 
the soft breaking parameters. The minimal realization, known as the Minimal 
Supersymmetric Standard Model (MSSM), has been the subject of numer- 
ous theoretical studies, and also experimental searches, for over two decades. 
While no experimental evidence for SUSY exists, the MSSM retains a pre- 
eminent status among models of TeV-scale physics in part through several 
indirect virtues, e.g. gauge coupling unification and a "natural" dark matter 
candidate. For full details of the MSSM spectrum and the parametrization of 
the soft-breaking terms, we address the reader to any of the comprehensive 
reviews on MSSM phenomenology [94]. 

The unbroken sector of the MSSM contains, besides the gauge interactions, the 
Yukawa couplings parametrized by 3x3 Yukawa matrices in flavour space, Y u , 
Y d and Y e . These matrices source the tree- level masses of matter fermions, 



where (Hi) and (H 2 ) are two Higgs vacuum expectation values related to the 
SM Higgs v.e.v. via (H2) 2 + (Hi) 2 = v 2 /2. In SUSY models, anomaly cancel- 
lation in the Higgsino sector requires the introduction of at least two Higgs 
superfields as above. In addition to Yukawa couplings, the supersymmetry- 
preserving sector contains the so-called //-term that provides a Dirac mass to 
the higgsinos (the superpartners of the Higgs bosons) and contributes to the 
mass term of the Higgs potential, 



Vkiggs = m\\Hi\ 2 + m\\H 2 \ 2 + m\ 2 HiH 2 + |//| 2 (|#i| 2 + \H 2 \ 2 ) + ■■■ (4.82) 



where the ellipses denotes quartic terms fixed by supersymmetry and gauge 
invariance [94]. to 2 , to 2 , and m\ 2 are soft-breaking parameters that may attain 




(4.81) 
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negative values thus driving electroweak symmetry breaking. By suitable phase 
redefinitions of Hi and H 2 , one can restrict to real Higgs v.e.v.s and introduce 
the parameter, tan/3 = (H 2 )/(H 1 ). 

Among the remaining soft-breaking parameters one has gaugino mass terms 
and squark and slepton masses, 

1 i=l,2,3 S=Q,U,D,L,E 

where Aj are the gaugino (Majorana) spinors, with i labelling the correspond- 
ing gauge group, U(l), SU(2) or SU(3). Each gaugino mass Mj can be complex. 
The second sum spans all the squarks and sleptons and contains five Hermi- 
tian 3x3 mass matrices in flavour space. Finally, the soft-breaking terms also 
include three-boson couplings allowed by gauge invariance, such as QH 2 A U U, 
that are called A-terms and are parametrized by three arbitrary complex ma- 
trices A u , Ad and A e . In the construction above we have limited the dis- 
cussion to the P-parity conserving case, which only allows an even number 
of superpartners in each physical vertex, and is imposed to reduce problems 
with baryon number violation. Even with this restriction, if we count all the 
free parameters in this model we find a huge number, of 0(100), with a few 
dozen new CP-violating phases! Truncating this number is fully justified only 
within the context of a fully specified supersymmetry breaking mechanism, 
which may then enforce additional symmetries and relations among parame- 
ters. 

Without going into the details of the dynamics behind SUSY breaking, it 
will be enough for our purposes to simply assume that the following, very 
restrictive, conditions are fulfilled: 

M| = m|l; for S = Q,U, D, L, E, "degeneracy" 

Aj = AiYi] for i = u, d, e, "proportionality", (4.84) 

Strictly speaking, such conditions can only be imposed at a specific normal- 
ization point above the weak scale, as the renormalization group evolution of 
the MSSM parameters will modify these relations. Moreover, these conditions 
can only be imposed with limited precision at this scale due to threshold ef- 
fects. Nonetheless, such a restrictive flavour universality ansatz in the scalar 
mass sector, and proportionality of the trilinear soft breaking terms to the 
Yukawa matrices has the utility that it greatly reduces the number of indepen- 
dent soft- breaking parameters. Even so, a significant number of CP- violating 
phases remain, e.g. 

Arg^M.m^), Arg(A / M*), Arg (MjM*), Aig(A f A* fl ). (4.85) 
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Fig. 7. One-loop SUSY threshold corrections in the down quark sector induced by 
a gluino-squark loop. On the left, a threshold correction generating Im(md), while 
on the right the analogous diagram for the EDM. The CP-violating source enters 
via the highlighted vertex, squark-mixing in the present case. 

Going to an even more restrictive framework, by assuming a common phase 
for the gaugino masses and another common phase for A$ reduces the number 
of independent CP violating parameters to two. Using phase redefinitions, one 
can choose the phase of the gaugino mass to be zero, and use 9a = Arg(A) 
and Op = Arg(/i) as the basis for parametrizing CP violation. 

It has been known for over twenty years that even in the absence of new 
flavour physics, large EDMs can be induced at the one-loop level within a 
single genaration [95,96]. Indeed, one would anticipate large EDMs as both of 
the reasons that rendered cL^Skm) very small, namely high loop order and also 
mixing angle/Yukawa coupling suppression, are not present for EDMs induced 
by the phases of the soft-breaking parameters. 

Figure 7 exhibits examples of one-loop diagrams at the supersymmetric thresh- 
old that generate non-zero contributions to the CP-odd Lagrangians (2.14) 
and (2.16). If we leave aside the problematic s-quark CEDM, then at one- 
loop we can concentrate on diagrams involving just the first generation of 
quarks and leptons. Within the parametrization described above, the phases 
residing in /1 and A permeate the squark, selectron, chargino and neutralino 
spectrum, which in the mass eigenstate basis translates into complex phases in 
the quark-squark-gluino and fermion-sfermion- chargino (neutralino) vertices. 
To make this explicit, for a moment let us truncate the flavour space to one 
generation and write down the expression for the 2x2 <i-squark mass matrix 
at the electroweak scale in the basis of di and cIr, 



Mi 



I 

V 



m 2 Q + 0{v 2 ) —md(p tan (3 + A* d) 
-m d ([i* tan/5 + A d ) m 2 D + 0(v 2 ) 



(4.86) 



where we further assume that the soft masses rriQ and m 2 D are large relative 
to the weak scale, and thus we can ignore subleading 0{y 2 ) corrections to 
the diagonal entries. Similar expressions can be written for the selectron mass 
matrix with the obvious substitutions in (4.86), and for the u squark, where 
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in addition one has to exchange tan (3 by cot (3. In the generic case of three 
generations, M 2 becomes a 6x6 matrix with 3x3 blocks which are tradition- 
ally called M| L , M\ R , M RL and M RR . For our purposes, the crucial terms in 
(4.86) are the off-diagonal components, 



(Mj) LR = -m d (fi tan (5 + A* d ) (4.87) 

which contain the CP-odd phases. By virtue of being proportional to the 
small mass m d , such a term can be treated as a perturbation and accounted 
for by an explicit mass insertion on the squark line, as in Figure 7. Note that 
the natural range for tan/5, in the interval between 1 and 60, allows for a 
significant enhancement of the /i-dependent term in (4.86). The phase of /i 
also modifies the spectrum of charginos and neutralinos, the mass eigenstates 
of the superpartners of the Z, W, 7 and Higgs bosons. 

With this notation in hand we see that, for example, the squark-gluino loop 
diagram generates an imaginary rf-quark mass correction that contributes to 



a s M 3 (n tan (3sm 6» M - A d sm 6 A ) Tfnjf , 

Im m d = -m d - — J(M 3 , m Q , m D ). (4. 

3vr M q 

The loop function I is normalized in such a way that I(m,m,m) = 1; its 
exact form (see e.g. [81,97]) is not important for our discussion. The ratio 
Im (m d )/m d , along with contributions from other quark flavours, represent a 
one-loop renormalization of the #-term. It is important to observe that it only 
depends on the SUSY mass ratio and thus does not decouple if A, fi, M3, tuq^d) 
are pushed far above the electroweak scale. Applying the bound on the 6*-term 
to the combined tree level and one-loop results (4.88), with degenerate SUSY 
mass parameters as above, we find |# t ree + 10 _2 5cp| < 10~ 9 , where #cp is a 
linear combination of sin# M and sin 6?^ with 0(1) coefficients. If there is no 
axion and 9 tTee vanishes instead by symmetry arguments, it follows that the 
phases of the soft-breaking parameters must be tuned to within a factor of 1CT 7 
in order to satisfy the EDM bounds. Therefore, an incredibly tight constraint 
on the phases of the SUSY soft-breaking parameters can be obtained in models 
which invoke high-scale symmetries to resolve the strong CP problem. 

However, if the PQ symmetry removes the #-term, such radiative corrections 
to 9 have no physical consequences, and the residual EDMs are determined 
by higher dimensional operators. The relevant expressions for the one-loop- 
induced di and di contributions can be found in e.g. [97]. Here we would just 
like to demonstrate the main point implied by these SUSY EDM calculations 
in a simplistic model in which all soft-breaking parameters are taken equal 
to a unique scale MgusY at the electroweak scale, i.e. Mj = mq = m D = 
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• • • — \fx\ — \Ai\ = Msusy- Working at leading order in f 2 /M| USY , we can 
then present the following compact results for all dimension 5 operators (with 
q = d,u), 

-t = M(' sin e (i [ tan/ j]±i - si n9 4 ) +0(gl,gl), (4.89) 

^ = ^ ( sin ^[tan/3]* 1 - sin 9 A ) + 0(gl gl). 

The notation [tan/3] 1111 implies that one uses the plus(minus) sign for d(u) 
quarks, gi are the gauge couplings, and e u = 2e/3, = — e/3. For the quarks 
we quoted the explicit result only for the gluino-squark diagram that dom- 
inates in this limit. All these contributions to di are proportional to a 
universal combination corresponding to the generic dipole size, 



« = 9 m \ = 1.3 x l(T 25 cm x _^L_ (i^V (4.90) 
16vr 2 M s 2 USY lMeV VM SU sy/ 

which varies by a factor of a few for i = e,d,u depending on the value of the 
fermion mass. The perturbative nature of the MSSM provides a loop suppres- 
sion factor in (4.90) so that Kj is about two orders of magnitude smaller than 
the estimate (4.80). Correspondingly, the reach of the current EDM constraints 
in SUSY models cannot exceed the scale of a few TeV. 

In (4.90) the quark masses should be normalized at the high scale, Msusy- To 
make the explicit connection with the dipole operators in (2.16), the results of 
Eq. (4.89) should be evolved down to the low-energy normalization point of 1 
GeV using the relevant anomalous dimensions (see e.g. [81]). Plugging these 
results into the expressions for d n , g?ti and o?Hg and comparing them to the 
current experimental bounds, we arrive at a set of constraints on 9a and 9^ 
depending on M S usy and tan/3. In Figure 8, we plot these constraints in the 
{Qufi a) -plane for M S usy = 500 GeV and tan/3 = 3. The region allowed by the 
EDM constraints is at the intersection of all three bands around 9 a = 9^ = 0. 
One can observe that the combination of all three constraints strengthens the 
bounds on the phases, and protects against the accidental cancellation of large 
phases that can occur within one particular observable. The uncertainty in the 
QCD calculations of g^NN and the nuclear calculation of S(g^ N ) discussed 
earlier may affect the width of the du g constraint band, but do not change its 
slope on the (6^, 9a) plane. 

Before we review the most common approaches to address the "overproduc- 
tion" of EDMs in supersymmetric models, for completeness, we will briefly 
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Fig. 8. The combination of the three most sensitive EDM constraints, d n , dxi and 
dug, for Mstjsy = 500 GeV, and tan/3 = 3. The region allowed by EDM constraints 
is at the intersection of all three bands around 9a = = 0. 

discuss some of the additional contributions which become important when 
tan j3 is large, a regime favoured for consistency of the MSSM Higgs sector 
with the final LEP results [98]. One simple observation is that the EDMs 
of down quarks and electrons, induced by 9^ at one-loop, grow linearly with 
tan/3 (4.89). However, at the two- loop level, there are additional contributions 
from the phase of the A-parameter which may also be tan /3-enhanced [99] . 
A typical representative of the two-loop family is presented in Figure 9. At 
large tan (3 the additional loop factor can be overcome, and these two-loop 
effects have to be taken into account alongside the one-loop contributions in 
(4.89). For example, the stop-loop contribution to the electron EDM in the 
same limit of a large universal SUSY mass is given by 



sin(^ + ^)tan/3, (4.91) 

where uia is the mass of the pseudoscalar Higgs boson, that we took to be 
smaller than Mstjsy, Y t — 1 is the top quark Yukawa coupling in the SM, and 
K e ~ 0.6 x 10 _25 cm. For very large values of tan/3 additional contributions 
from sbottom and stau loops, which are enhanced by higher powers of tan/3, 
also have to be taken into account [81,99]. 

Finally, the second, and in some sense more profound change is that at large 
tan/3, the observable EDMs of neutrons and heavy atoms receive significant 
contributions not only from the EDMs of the constituent particles, e.g. d e and 
d q , but also from CP-odd four-fermion operators [100]. The relevant Higgs- 
exchange diagram is shown in Figure 9. The CP violation in the Higgs-fermion 
vertex originates from the CP-odd correction to the fermion mass operator in 
Figure 7. These diagrams, since they are induced by Higgs exchange, receive 
an even more significant enhancement by (tan/3) 3 . In the same approxima- 
tion as before, the value of the thallium EDM induced by this Higgs-exhange 
mechanism, and normalized to the current experimental limit, is given by 
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Fig. 9. Additional corrections to the EDMs. On the left two-loop Barr-Zee type 
graphs mediated by a stop-loop and a pseudoscalar Higgs, while on the right 
we have a Higgs- mediated electron-quark interaction Cde with CP violation at 
the Higgs-quark vertex. There is a second diagram with CP-violation at the 
Higgs-electron vertex mediated by H . 

dm tan 3 /3 /100GeV\ 2 



[dm] cxp 

330 V mi / 



sin ^ + 0.04 sin(0 u + y 



(4.92) 



Notice that this result does not scale to zero as Mstjsy — * oo. Although just 
an O(10 -3 — 10~ 2 ) correction for tan/5 ~ 0(1), these Higgs-exchange contri- 
butions become very large for tan/? ~ O(50) [100,101,81] (see also [102]). 



4.2.1 The SUSY CP problem 

Figure 8 exemplifies the so-called SUSY CP problem: either the CP- violating 
phases are small, or the scale of the soft-breaking masses is significantly larger 
than ITeV, or schematically, 



\IVlf 



SUSY 



The need to provide a plausible explanation to the SUSY CP problem has 
spawned a sizable literature, and the following modifications to the SUSY 
spectrum have been discussed. 



• Heavy superpartners. If the masses of the supersymmetric partners exhibit 
certain hierarchy patterns the SUSY CP problem can be alleviated. One 
of the more actively discussed possibilities is an inverted hierarchy among 
the slepton and squark masses, i.e. with the squarks of the first two genera- 
tions being much heavier than the stops, sbottoms and staus, ie. (M|)y ^> 
(M|) i3 , (M|) 33 , where i,j = 1,2 is the generation index [103]. It is prefer- 
able to have masses of the third generation sfermions under the TeV scale 
because they enter into radiative corrections to the Higgs potential, and 
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making them too heavy would re-introduce the fine-tuning of the Higgs 
mass whose resolution was one of the primary motivations for weak-scale 
SUSY. Such a framework suppresses the one-loop EDMs which become im- 
measurably small if the scale of the u and d squarks is pushed all the way to 
~ 50TeV, as suggested by the absence of SUSY contributions in Am K ( B )- 
This does not mean, however, that the EDMs in such models become com- 
parable to ^(^km)- Indeed the two-loop contributions to di and w involving 
the third generation sfermions are not small in this framework, and indeed 
are at (or sometimes above) the level of current experimental sensitivity. 
Also, this means of suppressing the EDMs would not necessarily work in 
the large tan/3 regime where Higgs exchange may induce a large value for 
Cs that is not as sensitive to Msusy as the EDM operators. We note that 
future improvements in experimental precision will allow a stringent probe 
of such scenarios. 

• Small phases. A rather obvious possibility for suppressing EDMs is the as- 
sumption of an exact (or approximate) CP symmetry of the soft-breaking 
sector. This is essentially a "model-building" option, and various ways of 
avoiding the SUSY flavour and CP problem in this way have been suggested 
in the past fifteen years [104-106] . The idea of using low-energy supersym- 
metry breaking looks especially appealing, as it can also help in constructing 
an axion-less solution to the strong CP problem [48]. If the CP-odd phases 
in the soft-breaking sector are exactly zero and the conditions (4.84) are 
imposed exactly at the unification scale as a constraint on the high scale 
model, we can ask about the scale of EDMs induced by SUSY diagrams 
due purely to S KM - Since such an MSSM framework would possess the same 
flavour properties as the SM, one expects proportionality to the same CP- 
odd invariant combination of mixing angles, namely Jcp, and suppression 
by differences of Yukawa couplings [96] . Then it is easy to understand that 
the superpartner contributions to the down quark (C)EDM will necessarily 
be suppressed by the equivalent of 5 C p ~ JY^ ~ 10~ 9 , which is again six 
to seven orders of magnitude below current experimental capabilities, and 
thus not significantly larger than the EDMs induced in the SM [107]. 

• Accidental cancellations. Another possibility entertained in recent years 
[97,108] is the partial or complete cancellation between the contributions 
of several CP-odd sources to physical observables, thus allowing for 5 C p ~ 
0(1) with M SUSY < 1 TeV. Since the number of potential CP-odd phases is 
large, and the superpartner mass spectrum is clearly unknown, one cannot 
exclude this possibility in principle. However, as we illustrated in Figure 
8, d n , d-r\ and rf Hg depend on different combinations of phases, and the 
possibility of such a cancellation looks improbable. A more thorough explo- 
ration of the MSSM parameter space in search of acceptable solutions that 
pass the EDM constraints was performed in [109,110], and in the absence 
of additional parameter tuning did not identify any significant regions of 
cancellation. 
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• No electroweak scale super symmetry. Of course, there is always the possi- 
bility that other mechanisms (or no easily identifiable mechanism at all) lie 
behind the gauge hierarchy problem and the SM is a good effective theory 
valid up to energy scales much larger than 1 TeV. In this case there is no 
SUSY CP problem by definition. One of the recently suggested scenarios 
[111] exploits the possibility of a large number of electroweak vacua to invoke 
anthropic reasoning for selecting the "right" vacuum, thus side-stepping nat- 
uralness arguments for expecting new physics at the weak scale. Ref. [Ill] 
assumes that all the scalar superpartners are very heavy, but leaves gauginos 
and Higgsinos under a TeV, in order to preserve gauge-coupling unification 
and a dark matter candidate. This eliminates the one-loop induced EDMs, 
but leaves room for two-loop contributions [99,112] generated by chargino 
loops via a diagram similar to that shown in Figure 9 with A replaced by the 
light Higgs. This scenario can also be probed with the predicted sensitivity 
of future EDM experiments. 



4.3 SUSY EDMs from flavour physics 



EDMs can also serve as a sensitive probe of non-minimal flavour physics. In- 
deed, the assumptions of proportionality and universality in the soft-breaking 
sector (4.84) at a given high-energy scale are highly idealized, and are not 
expected to hold with arbitrary precision. In this subsection, we would like to 
show that EDMs are sensitive to flavour-changing terms in the soft-breaking 
sector, and provide significant constraints on SUSY models with non-minimal 
flavour structure. 

For concreteness, let us assume that (4.84) holds approximately, and the per- 
turbations are small. Around the electroweak scale, and in a basis with diago- 
nal quark mass matrices, the soft-breaking mass matrices can be approximated 

as 



M| = diag(m| n , m| 22 , m| 33 ) + SM 2 Sij , (4.94) 

where, as before, S labels the different squarks and sleptons, and % ^ j. Using 
this approximation, we can calculate the contributions to the relevant observ- 
ables using <5Mj^ as a perturbastion via insertions along the squark line, as 
in Figure 10. 

Calculating the gluino one-loop diagram in the approximation of equal SUSY 
mass scales, (M|)jj = Mf = |/i| 2 = M| USY , we arrive at the following result 
for the d quark EDM, and the imaginary correction to the d quark mass, 
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Fig. 10. Contribution of flavour changing processes to the <i-quark EDM. The middle 
insertion on the sfermion line corresponds to LR mixing proportional to m&; the 
insertions on the left and on the right correspond to flavor transitions in the LL 
and RR squark mass sector. 

r j a s tan (5 
lmm d = -b 131 x m b — — 

l07T 
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4 = 4xe "^ (495) 

where 5f 31 denotes the following CP-odd dimensionless combination, 



lm(5M" Ql3 e^5Ml 31 ] 

iw SUSY 



(4.96) 



In (4.95), for simplicity, we neglected the contributions from the A parameters, 
and retained only the mixing coefficients between the first and the third gen- 
erations. There are two important points about Eq. (4.95) that we should em- 
phasize here: <5i 3 i can be non-zero even if 6* M = 0, and both In^m^) and dd are 
enhanced relative to (4.88) and (4.89) by the large ratio (mj/m^) ~ 10 3 , which 
can compensate the suppression associated with flavour violation. In the case 
of u quark operators, this enhancement factor is even larger, m t /m u ~ 10 5 . 

As we have seen in the previous subsection, renormalization of 9 ~Im(m q )/m q 
can be very large, capable of producing bounds on Sf 31 and 5" 31 at the 10~ 9 level 
or better unless 9 is removed via PQ symmetry. In the latter case, using (4.95) 
and similar results in the lepton sector, one obtains the following sensitivity of 
EDMs to the above combination of flavour changing transitions on electron, 
u, and d quark lines for M SU sy = ITeV, 



s e 131 ~ icr 4 - icr 3 ; <^ 31 ~ icr 6 - icr 5 ; sf 31 ~ icr 4 - icr 3 . (4.97) 

Thus, EDMs independently provide very stringent constraints on the combined 
sources of flavour- and CP-violation in the soft-breaking sector. These con- 
straints are complementary to those coming from K and B meson physics and 
searches for lepton flavour violation. Note also that these calculations need to 
be modified in the large tan (3 regime to include Higgs-mediated contributions 
[113,100] which may dominate over the one-loop results for di (4.95). 
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It is important to realize that the apparent enhancement of EDMs in (4.95) 
by the ratios of heavy to light quark and lepton masses occurred because 
of the presence of flavour- changing terms in both LL and RR sectors of the 
squark/slepton mass matrices. Indeed, to make this point transparent we can 
write 



5 d 131 = Arg[(5 d 13 ) LL (5 d 33 ) LR (5 d 31 ) RR }, (4.98) 

in terms of "mass insertions" [114] <5/- = M~. ,/m|, which, although the distinc- 
tion is not crucial here, are usually defined in a slightly different basis to the 
one we have been using, m 2 - denotes here the average sfermion mass-squared. 
The status of the LL and RR insertions is in general rather different, and 
particularly so within the MSSM where the latter are essentially absent. 

To see this in more detail, we recall that flavour-changing terms in the LL 
sector are natural, as they are induced by renormalization group evolution of 
the soft-breaking parameters (see e.g. [115]) even if one assumes the conditions 
(4.84) at the unification scale. Starting from the universal boundary conditions 
(4.84) for all scalar masses, equal to Wq, and A parameters at some high-energy 
scale A uv , one can obtain an expression for Mq at a lower energy scale A, 
which at one-loop is given by 



Ml = ml 1 - ^i+d! ,„ (Mk) ( Y t y. + YjY„) + ■ ■ ■ (4.99) 

The ellipsis denotes "flavour-blind" contributions and also higher-order terms. 
Depending on the particular model of SUSY breaking Auv can be anywhere 
between a few tens of TeV and the Planck scale. The presence of both up 
and down Yukawa matrices in (4.99) guarantees the appearence of flavour- 
changing contributions in the LL entries of the squark mass matrices. At the 
superpartner threshold, A = Msusy, the flavour changing terms in the down 
squark sector will evidently be 



3m 



where V is the CKM matrix, and Y t is the top quark Yukawa coupling. If 
the scale A uv is very high, i.e. comparable to the Planck or GUT scale, the 
logarithm is large and can entirely offset the loop factor. Therefore, the natural 
size of the 13 entry in the down squark LL sector is ~ M| USY 14j — 0.01M| USY . 

The situation in the RR sector is completely different. There the absence of 
any Y u -dependence in the RG equations for M|, forbids the generation of sub- 
stantial flavour- changing transitions, unless the MSSM spectrum is modified 
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above certain energies so that the RG equations for the right-handed squark 
masses aquire flavour dependence. 

A number of SUSY scenarios have been proposed which describe plausible 
patterns of small deviations from (4.84), allowing for significant RR contribu- 
tions. In models with SO(10) unification and Ajjy > Aqut the running of the 
soft-breaking parameters extends above the unification scale, where the RG 
equations are modified by the presence of new field degrees of freedom. For 
SO(10) GUTs this modification introduces significant flavour dependence in 
the RR sector of squark and slepton mass matrices [116], even if the restric- 
tions (4.84) are imposed at the Planck scale. The resulting flavour-changing 
terms for down squarks are of the same order of magnitude as in the 

LL sector, leading to the prediction 8f 31 ~ 10~ 4 , which is right at the border- 
line of current experimental sensitivity (4.97), [117,118]. Similar effects may 
be generated by heavy sterile neutrinos. The light neutrino mass scale might, 
via the seesaw mechanism, be pointing to the existence of a new energy scale, 
Mr ~ YyV 2 j (0.1 — 0.001 eV) related to heavy sterile (or "right-handed") neu- 
trinos. If Kjjv is larger than Mr, the RG equations for sleptons will be modified 
above Mr with an effect similar to that above, namely a non-trivial flavour 
dependence will be imprinted on the slepton mass matrices. The importance 
of such an effect will depend on the size of the neutrino Yukawa couplings Y v 
and, with certain Yukawa patterns, an observable or nearly-observable elec- 
tron EDM might be induced [119]. Of course, if the scale of SUSY breaking 
is lower than Mr (or A GUT ) there are no significant consequences for EDMs 
unless one allows for other "diagonal" phases in this sector. 



5 Conclusions and Future directions 

Recent years have seen significant progress in the experimental tests of CP- 
violation in the Standard Model. Experimental verification of direct CP vio- 
lation in Kaon decay, and in particular the spectacular measurements of CP 
asymmetries for neutral B meson decays at BaBar and Belle have provided 
solid confirmation of the overall validity of the Kobayashi-Maskawa mech- 
anism. The current status of CP violation in flavour changing processes is 
such that (within errors) it does not necessitate the introduction of any addi- 
tional CP-violating sources. At the same time, there is ample (experimental) 
room for the existence of new CP- violating physics to which the K and B 
meson data is not sensitive. This concerns, primarily, CP violation in flavour- 
conserving channels. The existence of such new sources is hinted at, albeit 
indirectly, by the baryon asymmetry in the Universe. The search for CP vio- 
lation in flavour-conserving channels, and the search for EDMs in particular, 
should thus remain high on the priority list for particle physics. The strong 
suppression of EDMs induced purely by the Kobayashi-Maskawa phase, com- 
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bined with prospects for improving the experimental sensitivity, places EDM 
searches at the forefront in probing CP- violating physics beyond the Standard 
Model. 

Moreover, beyond their direct sensitivity, current (and future) null results 
for EDM searches also provide very powerful constraints on models for new 
physics. Indeed, as we have discussed, the sensitivity for example to CP vi- 
olation in the soft-breaking sector of SUSY models, allows us to probe soft- 
breaking masses as large as a few TeV. In this indirect sense, EDMs are often 
sensitive to energy scales beyond the reach of future collider experiments, and 
play a central role in the full suite of precision tests of the Standard Model. 
The scales probed by EDMs and by the constraints on flavour changing neu- 
tral currents are not too dis-similar, and the gap may continue to narrow with 
future progress in EDM searches. This only heightens the tension between the 
observed CP- violation in the flavour- changing sector and the lack thereof in 
flavour-diagonal channels, of which the strong CP problem is the most man- 
ifest example. We seem compelled to question whether CP and flavour are 
as intrinisically linked in general as they are within the Kobayashi-Maskawa 
model? This is one aspect of what one might hope would be answered by a 
general "theory of flavour". EDMs will clearly continue to provide a crucial 
probe in tackling this question. 

In this concluding section we would like to emphasize some directions on the 
experimental and theoretical front that are likely to bring future progress in es- 
tablishing the nature of flavourless CP-violation at and above the electroweak 
scale. 

• Experimental Developments 

There are a number of developments in experimental techniques to search for 
EDMs which promise to narrow the gap between the current limits and the 
KM background in all of the EDM classes discussed in this review. In par- 
ticular, novel techniques for storing neutrons in liquid helium, in progress at 
LANSCE in Los Alamos [120] and under investigation for an experiment at 
PSI, will help to improve the measurement of many fundamantal parameters 
in neutron physics, including d n . On another front, the suggestion that CP- 
violating effects can be significantly enhanced in exotic nuclei possessing an 
octopole moment [121] drives several experiments searching for EDMs of iso- 
topes of Ra and Rn [122] . We note that such experiments should be pursued 
primarily because of their potential to discover CP violation, as null results 
will not be as constraining as those from d n or du g due to large uncertain- 
ties in the calculations of nuclear matrix elements. In future measurements 
of paramagnetic EDMs, the resulting sensitivity to d e and Cs will be signifi- 
cantly improved via the use of paramagnetic molecules, such as YbF and PbO, 
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that can be polarized and thus allow a huge enhancement of the applied field 
[15,23]. The anticipated precision will allow probes of the electron EDM down 
to lCT 30 e cm or below [15,23]. In this regard, we should mention an interest- 
ing alternative approach, which involves the measurement of the tiny electron 
EDM-induced magnetic flux when an electric field is applied to a particular 
garnet crystal. This project is also already in development at LANSCE [123]. 

Perhaps the most interesting proposal of recent years is the new approach to 
measure the EDMs of charged particles in storage rings. An initial proposal 
to measure the muon EDM at the 10~ 24 e cm level [124] has since evolved 
into the idea to measure the EDMs of ions and nuclei [125], and in particular 
the deuteron EDM at the 10 _27 e cm level. Although clearly not a diamagnetic 
system, the deuteron EDM can be placed in the same category as it is primarily 
sensitive to the same nuclear scale CP-violating source, namely the CP-odd 
pion-nucleon couplings g n NN- More precisely, one finds [126-128], 



d D = (d n + d p ) - (1.3 ± 0.3 GeV-^e^, (5.101) 

where the second term generically dominates, leading to a similar depen- 
dence on CP-odd sources as for mercury. However, in comparison to mercury, 
the deuteron EDM proposal has the significant advantage of requiring only 
straightforward nuclear calculations, since it is a weakly bound state, which 
leads to the relatively precise dependence on g^ N in (5.101). Using the result 
(3.71) for 

9ttnni an d accounting for the subleading corrections, one obtains 

[127] 



d D ~ Q + _fe{d d - d d ) + 0(d d + d u , d u , d d ). (5.102) 

Thus a measurement of do at the ~ 10~ 27 e cm level would correspond to 
a sensitivity to light quark CEDMs at the level of 10~ 28 cm [127], which is 
at least two orders of magnitude better than the current limits from dn g . 
Above all else, these proposals to measure EDMs in storange rings deserve 
special considerations as they depart from the dominant philosophy that EDM 
measurements demand the study of neutral objects in parallel electric and 
magnetic fields. 

To place this activity in context, we should bear in mind that probably the 
most important single question for particle physics - the origin of electroweak 
symmetry breaking - will be subjected to serious experimental scrutiny with 
the Large Hadron Collider coming on line within a few years. Besides the dis- 
covery of the Higgs boson(s), it may provide an answer to the gauge hierarchy 
problem, and indeed uncover a plethora of new particles or resonances above 
the electroweak scale. EDM experiments, which might of course discover new 
physics before the LHC switches on, can subsequently play a complementary 
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role in providing constraints on (or signatures of) CP-violating couplings (e.g. 
in the Higgs sector of the MSSM). The projected level of sensitivity in coming 
years will be more than competetive in this regard with collider probes. More- 
over, strangely enough, the absence of new physics (beyond the Higgs - or 
whatever might play this role) at the TeV scale would not remove motivations 
for EDM searches. Indeed, as we argued in this review, EDMs are sensitive to 
CP violation at multi-TeV scales, and thus represent one of the few classes 
of low-energy precision measurements that are sensitive to such high-energy 
scales. 

Another important experimental direction relevant to CP- violating physics is 
the search for axions. As we reviewed, one of the more natural resolutions of the 
strong CP problem predicts the existence of a light pseudo-scalar particle, the 
axion. The developments of recent years in cosmology have lent considerable 
weight to the presence of a non-baryonic cold dark matter component of the 
energy density in the Universe. Although the popularity of supersymmetric 
models continues to focus attention on the lightest supersymmetric particle 
(or LSP) as a natural dark matter candidate at the weak scale, axions with 
a coupling f' 1 below its astrophysical bound in fact still represent a viable 
alternative, thus providing additional motivation for the continuation of axion 
searches. 

• Theoretical Developments 

On the theoretical side, beyond questions of the precise generation mecha- 
nisms of CP-odd sources in specific new physics models, it is clear that the 
primary limitation on the full application of the observational bounds arises 
through the limited precision of QCD and nuclear calculations. Perhaps the 
most afflicted quantity at present is the CP-odd pion-nucleon constant, as in- 
duced in particular by the CEDMs of light quarks. As we have discussed, this 
is a fundamental parameter controlling the level of the constraints imposed by 
diamagnetic atoms, which can currently be calculated only to limited precision 
due to large cancellations in the relevant nucleon matrix elements. Another 
important issue concerns the strange quark CEDM contribution both to g^NN 
and the neutron EDM d n , and whether or not it is underestimated in the 
leading-order sum-rules analysis [59,129]. It would clearly be worthwhile to 
revisit these aspects. However, it seems likely that significant quantitative 
progress will come only from ab initio lattice calculations. This is a very chal- 
lenging task, since a viable lattice calculation would necessarily have to respect 
chiral symmetry both at the level of quarks and gluons and also among the 
observable matrix elements between the hadronic states, since this is the un- 
derlying reason for the suppression of d n (6) by m* and the partial suppression 
of d n (d q ). To that end, it will be important to implement a calculation dis- 
playing all the required symmetries, and in this sense d n {6) would be a good 
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place to start, as many features of the answer, such as the dependence on 
m, and on 6 = 9 + arg detM q , are enforced by symmetry allowing for inde- 
pendent checks of the calculation. On the nuclear side, we noted that recent 
reanalyses of the Schiff moment indicate that various many-body effects, e.g. 
polarization, can be significant and thus further progress in this area would 
assist significantly in improving the quality of constraints on g w NN in differ- 
ent isospin channels. It will also be important, in guiding future experimental 
ideas, to clarify the size of the enhancement of CP violation in exotic nuclei 
with octupole deformations. 

In conclusion, the limits on flavour- diagonal CP- violation produced by the null 
results of existing EDM searches already provide strong constraints on new 
physics at and above the electroweak scale. Developments in coming years 
promise to provide us with a wealth of new information about the nature of 
CP violation and TeV-scale physics, complementary to studies of electroweak 
symmetry breaking at colliders and flavour studies with K and B mesons. 
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